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Measuring  Non-Stat ionar i t y in  the  Parameters 
of  a Linear  Model  witli  Applications 
to  Asset  Returns 


by 

Melvin  J . Hinich  and  Richard  Roll 
\ ABSTRACT 

Tlie  purpose  of  this  paper  is  to  introduce  a technique  for  measuring 
non-stationarity  in  linear  models.  The  technique  is  conceptually  simple. 

It  can  be  adapted  easily  to  a wide  variety  of  problems  and  improved 
upon  in  obvious  ways  (whicli  we  will  mention).  Its  restriction  to  linearity 
is  probably  not  important  to  the  practitioner  because  the  overwhelming 
majority  of  commonly  used  techniques  are  linear;  Lxamples;  the  sample 
mean,  regression,  analysis  of  variance.  method  produces  a time 

path  of  a linear  model's  coefficients  and  it  provides  the  capability  to 
assess  the  statistical  significance  of  the  non-stationarity.  Fur  t liermi're , 
this  is  done  under  very  weak  assumptions  concerning  the  probability 
distribution  that  generated  the  data.  No  particular  generating  process 
need  be  assumed.  .Such  a robust  method  will  cleai  ly  not  be  optimal  lor 
every  application.  An  investigator  who  knows  the  generating  process, 
or  is  willing  to  assume  that  he  knows  it,  will  be  able  to  find  a 
specialized  technique  better  adapted  to  his  case;  but  for  those  <*t  us 
with  a lower  degree  of  skill  or  of  arrogance,  the  method  to  be  described 
here  has  much  practical  value. 


i 


I.  INTRODL'CnnN 


Stationarlty  over  tine  Is  aa  Important  unsolved  problem  In  the  application 
of  mathematical  models.  In  banking  and  finance,  the  problem  is  critical,  be- 
cause statistical  techniques  are  frequently  applied  to  time  series  that  extend 
backward  over  many  periods.  Virtually  all  analyses  of  time  series,  whether 
they  involve  the  simple  calculation  of  a mean  (such  as  a bank  s average  daily 
foreign  currency  position)  or  very  complex  calculations  (such  as  forecasts  of 
Inflation  generated  by  a 300  equation  econometric  model)  are  subject  to  the  non- 
atatlonarlty  problem.  There  is  usually  no  theoretical  reason  nor  practical 
guarantee  that  the  data  used  for  estimation  were  generated  by  a stationary  pro- 
cess. 

The  purpose  of  this  paper  is  to  introduce  a technique  for  measuring  non- 
statlonarlty  in  linear  models.  The  technique  is  conceptually  simple.  It  can 
be  adapted  easily  to  a wide  variety  of  problems  and  improved  upon  in  obvious 
ways  (which  we  will  mention).  Its  restriction  to  linearity  is  probably  not 
Important  to  the  practitioner  because  the  overwhelming  majority  of  commonly- 
used  techniques  are  linear;  Examples: the  sample  mean,  regression,  analysis  of 
variance. 

Our  method  produces  a time  path  of  a linear  model's  coefficients  and  it 
provides  the  capability  to  assess  the  statistical  significance  of  the  nonsta- 
tionarity.  Furthermore,  this  is  done  under  very  weak  assumptions  concerning 
the  probability- distribution  that  generated  the  data.  No  particular  generating 
process  need  be  assumed.  Sucha  robust  method  will  clearly  not  be  optimal  for 
every  application.  An  investigator  who  knows  the  generating  process,  or  is  willing 
to  assume  that  he  knows  it,  will  be  able  to  find  a specialised  technique  better 
adapted  to  his  case;  but  for  those  of  us  with  a lower  degree  of  skill  or  of 
arrogance,  the  method  to  be  described  here  has  much  practical  value. 

W*  will  Illustrate  the  method  by  using  a well-known  linear  model  of  asset 
returns,  the  Index  equation. 
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R»”Qf+/}R+C  ( 

J.t  J.t  m.t  J.t 

In  which  R is  the  observed  rate  of  return  on  asset  J during  time 

period  t,  R is  the  rctunr  on  a market  index,  ff  is  a parameter 
m,t 

which  measures  the  risk  of  asset  j , Is  another  parameter  whlcli 

■'  1/ 

has  been  given  several  interpretat  ions  by  different  authors,  and 
Is  a stochastic  disturbance  term  with  zero  mean. 


Since  excellent  expositions  of  the  meaning  of  (1)  arc  availaf.lc  elsewhere, - 
we  will  not  give  yet  another  one.  Suffice  it  note  that  (1),  or  closely  related 
variants,  are  widely  applied  in  the  financial  industry.  In  portfolio  management , 
there  are  probably  no  current  paradigms  more  important  than  (1).  It  is  an  aid 
to  decision  making  in  risky  environments.  It  has  been  used  for  assessing  port- 
folio  performance,  for  measuring  the  Impact  of  various  events  on  security  prices, - 
for  measuring  rates  of  discount  to  use  in  the  valuation  of  uncertain  cash  flows, 
and  for  other  purposes. 

In  most  applications,  the  parameters  j and  ^ j are  assumed  constant. 
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TTJierc  Jrc  however,  three  possible  sources  of  non-ststlonsrlty  which 
could  have  s detriments!  Impact  on  the  quality  pf  estimates  derived  frcjro 
fitting  (1)  to  time  series  of  returns,  while  assuming  stationary  paramet- 
ers, First,  there  is  no  theoretical  reason  for  the  "risk"  parameter  j 

to  be  an  Intertemporal  constant.  In  fact,  this  parameter  is  widely- 

acknowledged  to  depend  on  factors  that  are  known  to  change;  (such  as 

4/  1 \ 

the  debt-equity  ratio  of  the  firm—  when  asset  J is  a common  stock), 
and  the  absolute  risk  aversion  of  Investors  (which  is  likely  to  change 
with  the  aggregate  level  of  wealth) 

Another  source  of  non- stationarity  is  in  the  distribution  of 
the  disturbance  term  Intended  to  measure  all  influences  on  the  return  of 
asset  J that  are  unrelated  to  the  Index,  m.  — ^ There  exista  no  satis- 
factory theory  about  the  distribution  of  dj  aside  from  the  simple  assertion 
that  Its  mean  la  rero.  In  any  finite  sample,  of  course,  the  mean  of  6 ^ 
can  be  non-iero  In  so  far  as  dls-equllibrium  conditions  occur  temporarily. 
Also,  the  higher  moments  of  the  distribution  of  dj  need  not  be  constant. 

For  example,  the  variance  of  dj  could  fluctuate  drastically  and  cause 
considerable  difficulty  in  econometric  estimation  of  (1). 

Finally,  non-statlonar Ity  can  be  Induced  In  the  model  and  in  estimates 
of  its  parameters  if  the  particular  form  (1)  is  mis-spec if led . This  Is  the 
situation  %«hen  the  true  state  of  nature  requires  explanatory  factors  In 
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addition  to  the  market  return, ■"  1^.  Since  some  recent,"" 
and  some  not  so  recent  evidence  argvres  for  the  Inclusion 

of  other  factors,  thrs  source  of  non-stat ionarity  should  not  be  ruled 
out  in  advance.  We  shall  not  cxainlnc  the  question  of  v;hlch  pro- 
posed factors  should  be  included.  Omitted  variables  cause  the  parameters  of 
(1)  to  become  functions  of  time.  This  could  occur  an>'way,  even  after  cotiec- 
ting  the  mis-spcci flcation , so  ve  think  little  is  lost  and  much  clarity  is 
gained  by  developing  techniques  for  measuring  non-stationarity  with  the 
model  in  its  simple  form. 


An  Outline  of  V.’hat  Follows 

In  the  two  sections  following,  we  will  discuss  techniques  for  avoiding 
some  of  the  pitfalls  caused  by  non-stationarfty.  llic  next  section  (11)  vill 
explain  the  technique  of  robust  regression;  that  is,  of  an  improved  re- 
gression method  for  a model  such  as  (1)  vthon  the  disturbances  do  not  have 
all  the  standard  spherical  Gaussian  properties  so  familiar  from  econometrics 
texts . 

There  is  no  doubt  that  asset  returns,  and  the  estimated  residuals  from 

regressions  with  asset  returns,  do  not  have  standard  properties.  Asset 

returns  are  not  stationary  Gaussian.—^  There  is  a controversy  about 

whether  the  disturbances  of  (1)  are  drawn  from  a stationary  non-Caiissi an 

12  * 

probability  di.stributlon,  such  as  a stable  distribution,—  or  whether  tlicy 

5 3/  , 

conform  better  to  a Gaussian  process  with  non-stat jonary  parameters; — but 
whatever  the  true  explanation  may  fbc  observed  distributions  strongly 

suggest  that  an  alternative  to  the  classic  least  squares  method  of  line 
fitting  may  give  better  results.  Wo  will  pro/ldc  evidence  that  "robust 
regression"  does  indeed  perform  better  than  ordinary  least  squares. 


l« 
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In  eection  llifcc,  we  descrlliC  ttio  second  part  of  n mclbodolocy  for  re- 
Crcsslon  with  n non-statlonary  linfar  model.  It  utilizes  orthogonal  poly- 
nomials ot  time  to  track  the  paths  tnson  by  the  coefficients,  tli.-  "risk" 
coefficient  , an  ! the  intercept  paraMeU  r , . Section  four  gives 

the  empirical  results  obtained  from  a sample  of  U.S.  asoets. 

In  summary,  the  basic  model  (1)  will  be  fit  to  data  while  or  , 

P and  the  distribution,  f (C  ),  of  the  disturbance  term  are  fit 

jit  j|ijib 

to  explicit  functions  of  time.  This  is  a very  robust  specification  and 
It  has  the  potential  to  side-step  many  of  the  troubling  theoretical  and 
econometric  problems  previously  mentioned.  For  example,  if  the  "true" 
state  of  nature  requires  another  factor  (such  as  in  the  models  of  Black 
[1972]  or  of  Merton  [1973]),  the  omission  of  this  variable  is  cconometrically 
corrected  in  our  specification  by  allowing  the  intercept  to  vary  in  time. 
Similarly,  if  the  simple  Sharpe-Llntncr  model,  hut  with  non-constant  para- 
meters, is  the  "true"  state  of  nature,  our  specification  will  work.  Ey 
using  model  (1)  we  are  not  obliged  to  take  sides  on  the  question  of  which 
particular  theory  Is  "true"  because  this  specification  will  be  an  ap- 
proximation to  all  of  the  currently  suggested  theories  (for  a particular 
data  sample) . 


6 


TJ.  ROMISr  RKf.UKS.SlO.-) 

The  parameters  of  a linear  model  are  often  cstlriatccl  by  the  rethod  of 
ordinary  least-squares.  (OI.S)  which  is  sensitive  to  larco  values  of  the 
additive  error  terms.  Various  alternative  nctliods  have  been  proposecl  for 
obtaining  regression  estimates  v;hich  are.  insensitive  to  large  disturbances 
and  have  known  sampling  properties,  at  least  asymptotically.  (See  bickcl 
(19731  and  Huber  [1973]).  The  Iterative  algorith;:;s  used  in  all  th.-se 
studies  have  the  disadvantage  that  they  require  a preliminary  reasonable 
estimate,  usually  OI.S.  Using  linear  prograirjning  to  minim- ?.c  the  sum  of 
absolute  errors  protects  agafnst  large  disturbances,  but  the  only  known 
sampling  results  were  found  by  artiffcial  data  studies  (Blattbcrg  and 
Sargent  (1971]).  Hinich  and  Talwar  [1975]  present  an  alternative  simple 
two  stage  procedure  and  study  Its  asymptotic  and  empirical  proponties  wlien 
the  disturbances  are  assvimed  to  be  independent  realizations  from  a sys.imettic 
stable  distribution. 

Let  us  present  the  Uinlch-Talwar  procedure  for  estimating  the  coeffic- 
ients a and  P in  the  simple  model. 
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been  eliown  to  be  "lobubL"  iUi;luibnnrrs  (Incliidlnf,  (-.ro-s  data 

errors)  provided  tlic  cocf  fi  clouts  fr  and /)  are  con‘.tani  over  time. 

Assuming  for  convenience  that  T Is  even,  divide  the  sample  into  1/2 
non-overlapping  groups  of  two  successive  observations  ^^t+J  ' 

V ).  t " 1 3 ....  T-1  . (Note  lli.at  in  nodels  with  k variable:.  fV.  >2), 

•'t+1  . . 1 

th«  sample  would  be  divided  into  T/k  non-overlapping  groups.) 

For  each  t,  compute  the  equation  of  the  line  connecting  the  pair  of  points; 

1 .e. , 

Vl^t  - ^ytdl 

« 

'^+1  - ^ 


and  0 


(and,  of  course,  there  would  be  a vector  of  k estic.atcd  coefficients  in  the 
k-varlable  case.) 

0 t •••!  0 ®te  independent  r andoui  va  r 1 ah  Ics  with  the  same  loca  1 1 on  , 0 , 
but  with  different  scales— ( (or  dispersion  parameters).  The  scale  of  0^  Is 


0^(0) 


.1^+1  - 


t - 1,3,...,  T-1 


Slmi larly , , . 

and  scales 


0^  are  independent  random  v.Tii.sblrs  with  lor.ition  a 


o^Ca) 


(ivird-.i) 


t » 3,3,..,,  T-1 


An  initial  estimate  of  0 is  given  by 


0 » Median  (0  . ,0  ^.j)  , 


1/2  -1 

which  l.s  a consistent  and  asympotical  ly  normal  estimator:  T 6.^,  (0-0) 


21 
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converge:;  In  Jlutrlbutlon  to  N(0,1)  as  T -►  ®,  where 


6^  ■=  4 ( 

T 2 ‘ 


rO/Y) 


I (p)  ]■ 


If  the  X^s  are  stochastic,  it  fullovs  from  (6)  that  for  large  T,  P is 
2 

approximately  ti(/l  fi  /T)  where 


, lA-1/2 

. I -1 

0 (K|x  - X I ) 

ro/Y) 


The  asymptotic  efficiency  of  p is  similar  to  that  for  the  median 
estimate  of  the  location  given  in  a random  sample  from  a parent  distri- 
bution. The  efficiency  of  the  first  stage  estimator  can  be  increased  by 
using  a truncated  mean,  (e.g.,Thc  mean  of  the  middle  25%  of  the  ordered  ' s , 
Instead  of  the  median  ^ • The  median  was  presented  here  because  its  asym- 
potic  properties  arc  easier  to  derive  and  express.  The  first  stage  estimator 
of  a has  similiar  properties. 

Once  the  first  stage  estimates  of  a and  /?  arc  computed,  the  residuals. 


■=  y,  - « -P\ 


are  ordered  and  used  to  compute  an  estimate  of  c , 


® 1.654  ^*^(.721)“  ®(.28T)^’ 

where  c.  T-.,,.  and  c,  at'e  the  order  statistic  estimators  of  the  28th 

(.til)  I . ^oi ; 

and  72nd  percentiles  of  the  distribution  of  c.  The  cstin\ator  s h.as  an 
asymptotic  bias  of  less  than  0.4%  for  .all  y ln"[l,2).  (See  Kama  .and  Soil 
[1968,  p/  823]). 


<*:w 
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The  6o;ni>Jc  le  llun  dcclmat  cd  hy  removing  nil  olir.c  rva  tloiis 
re'ipondlnr,  to  rci-ldiial;.  wliJch  art-  greater  In  absolute  value  than  nuir.o 
fixed  multiple  of  s.  Here,  we  used  ''is  as  tlic  cutoff.  This  value  was 
shoyu  in  the  monte  carlo  experiraents  of  Hinich  ond  Talvar  to  bring 
relatively  low  sampling  dispersion  to  the  resulting  coefficients  over  a 
wide  range  of  distributions  ot  the  raw  disturbances.  The  final  estimates 
of  a and  p arc  ttie  ordinary  leas t-squarcs  coefficients  computed  from  the 
remaining  observations.  They  arc  approxlr,.? tely  nomal  1 y-dintributed  since 
the  errors  in  the  renaining  observations  have  finite  variance.  As  we 
shall  see,  this  ability  to  "force''  the  second-stage  1 ca.s  I -squares  es  L iin.'it  < 
toward  normality  Is  of  crucial  importance  in  testing  for  their  non- st.'-tior,- 
arity.  Furthermore,  the  Hinich-Talwar  technique  will  provide  this  result 
for  any  atochastic  process  of  the  disturbances,  t,  not  just  for  stable 
non-Gaussian  disturbances,  A non-stationarity  in  the  dispersion  of  C,  for 
example,  will  bring  about  thic’xer  tails  in  the  sampling  distribution  of  C.  , 
By  using  this  robust  procedure,  however,  the  least-squares  estimates  com- 
puted under  these  conditions  will  al.so  be  forced  toward  normality. 


1 
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III . HorKi.i.i’.'o  ?:o,\-STATin:;AUTri’  ir:  ti'i-  i::'rrgr!  I'i . ':i'k  .si.npr.  onKFKir.i-t^jT 

A>iD  T!;i.  ' 

In  sejrcliliif;  for  ttic  bcHl  technlijuc  to  model  non-r. taLlonarl  ty  1 ii  tlio 
i 

I 

parpmcters  of  equation  (1),  there  arc  several  useful  Implications  of 
past  empirical  woi  1,  tliat  should  be  considered.  1 lie  previous  section  out- 
lined an  approach  for  accomodating  nen-st at lonar i t y In  the  dispersion  of 
the  disturbance  terra  . Here,  we  wish  to  present  a treatment  of  non- 

j 6 1 at  I onari ty  in  the  expected  value  of  C and  in  the  coefficients.  I 

I 

For  common  stocks,  several  studies  have  documented  temporary 
deviations  of  E(C)  from  zero.  This  seems  to  occur  as  a result  of  market 
dlscquilibria  when  new  and  unantl cl pated  information  is  received  by  in- 
dividual traders.  For  example,  the  mean  disturbance  term  is  significantly 
positive  in  the  w'eeUs  preceding  a stock  split— ^ In  another  example,  pcsi- 
' tlve  means  occur  before  an  announcement  of  spuriously  increased  earnings 

(caused  by  accounting  manipulations)  while  negative  means  follow  sub- 

1 fi  i' 

sequent  disclosure  of  the  spurious  nature  of  the  increase — '.  Non-zero 

mean  disturbances  have  been  associated  empirically  with  several  other  oc- 

19,  , . 20/  ^ , 

curences  such  as  secondary  offerings — and  dividend  increases — and  tneic 
are  undoubtedly  many  other  circumstances,  as  yet  not  documented,  which  have 

t 

21  / 

the  same  result. — 

d 

The  same  arguments  can  be  advanced  concerning  , the  intercept.  If  the  ^ 

Intercept  Is  interpreted  as  in  the  Sharpe-Lintner  model  as  a function  of  the 
riskless  rate  of  interest,  it  can  vary  intcrtemporally  too.  There  is  pro- 

1 

bably  not  as  much  solid  empirical  support  for  its  significant  variation;  but  | 

there  is  certainly  no  reason  to  assume  constancy  in  the  absence  of  evidence  | 

either  way.  In  several  otlicx  studies  (previously  mentioned),  the  intercept  ! 

• ■ 1 

has  been  presumed  to  vary  with  some  other  stochastic  return.  There  is  still 

I 
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a controversy  about  v!ilc!i  Ofisel  1)01011^’'^  In  the  model  and  since 

we  can  account  for  r.on- s t « 1 1 ona  r 1 1 y In  ^ by  an  altcinatlve  pro- 
cedure, we  perfer  to  not  measure  the  "lals.slny  fnrtot"  dirt-ctly.  sim- 
ple reflection  shows  that  (1)  permits  a methotl  for  naosuring  and 


correcting  non-st  at  ionnrl  ties  in  both  b(C^^)  and  in  the  iwlGbiiig  factor  '.or 
the  riskless  return)  at  the  same  time.  The  temporally  varying  Intercept 
can  simply  be  regarded  as  a mixture  of  all  these  influences. 

Vliich  time  dependent  functions  should  be  ctiosen  to  model  c/.,  and  p ? 

j » t J 1 1 

Several  alternatives  are  available:  One  of  the  techniques*  that  might  have 
been  selected  is  a niodi  f i cat  i on  of  the  Farley-lil  ni  t h (1970)  or  Fax  1 ey-Hi  ni  c.u- 
McGui  re  [1973]  procedures,  which  allow  the  coefficients  in  a rcg.ression 
equation  such  as  (1)  to  shift  once  during  the  observed  record.  Their 
technique  allows  the  (unknown)  date  of  the  shift  to  be  estimated.  Ihiu  woulT 
have  been  an  excellent  method  if  changes  in  the  coeffielents  of  (I)  ocruxred 
discretely  (and  only  once  dui'ing  the  pexiod  of  observation).  There  may  be 
cases  like  this.  In  fact,  such  a case  would  occur  for  a firm  thst  hud 
floated  a new  bond  issue  and  thereby  changed  crnr.i  dcrabl  y its  lever<>ge  ratio 
during  the  sample.  But  there  are  certainly  more  complex  possibilities  and 
prudence  required  a more  general  r.i'thoci. 

Another  alternative  is  seme  type  of  "random  coefficients"  tech- 
nique such  as  the  adaptive  regression  procedure  of  Cooley  and 


Prescott,  (1973).  In  their  method,  the  coefficients  arc  assumed  to  vary 
randomly  from  one  period  to  the  next  by  following  a random  xralk  from  their 
Initial  positions.  Again,  the  coefficients  of  sc'me  assets  may  bcliavc  this 
way  but  others  may  take  on  discrete,  determini' i ic  changes  or  fluctuate  In 
a predictable  pattern.  Besides,  the  methodology  described  belox'  is  able  to 
track  coefficients  which  actually  follow  a random  sequence.  Since  it  will 


i 

\i 
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track  other  sequences  too,  iiciic  as  aosi.  oy  xls  ooi;. 

Our  Idea  is  to  approximate  the  sequence  of  each  coefflcacnt  by  a func- 
tion of  time  whose  parameters  can  be  estimated  directly.  We  decided  to 


use  a function  of  Legendre  polynomials  for  this  purpose.  The  Legendre  poly- 
nomials are  only  one  of  a wide  varlery  of  fuurrlono  that  night  he  u'u'd  In  dif- 
ferent applications.  We  do  not  Imply,  hy  using  the  Legendre  polynomials  here, 
that  they  are  necessarily  Superior  to  Chebyshev  polynomials,  trigononet ric 
polynomials  or  many  others  that  could  be  used  to  approximate  any  arbitrary 
function  of  time.  However  the  Legendre  polynomials  are  easy  to  visualize 
and  they  provide  an  adequate  expositlonal  device  to  Illustrate  the  general  prin- 
cipals involved  In  our  technique.  * 

Given  any  polynomial  of  time,  the  time  paths  of  the  coe f f ic 1 ent s tx j ^ and 
can  approximate  as 


“j.r  Ji-o 


( - 1 ^ t 1)  (8) 


( - 1 ^ t ^ 1) 


(9) 


where  P^(t)  is  the  polynomial  of  order  i at  time  t and  j and  bjj  are  the  co- 
efficients of  the  ith-order  polynomial  for  the  intercept  and  slope  coefficients 
of  (1),  respectively.  In  this  case,  ? is  the  set  of  Legendre  polynomials  and 

the  first  five  of  these  arc  given  in  Table  1 below  and  illustrated  In  Figure  1. 
Units  of  time  are  chosen  so  that  -1  to  +1  spans  tlic  observed  number  of 


21/ 

natural  calendar  units — . This  particular  structure  is  used  because  the 

2a  / 

Legendre  polynomials  arc  mutually  orthogonal  on  the  interval  [-1,1 
Thus,  cacli  polynomial  can  be  introduced  as  a separate  variable  in  the  es- 
timating regression  without  having  to  worry  about  .mul ti-col li near! ty  ; (r’hich 
would  be  serious  problem  if  a non-or thogonal  polynomial  were  used). 


The  final  estimating  equation  is: 


OJm,t  ljm,t  2J  • m,t 
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ToMo  1 

Legendre  Polynomials  for  unit  wclgliL  and  range  (-1,  1). 

p^(t)  - 1 

Pj^(t)  ^ t 

P2(t)  " (3l^  - 1)/?. 

p^Ct)  “ (5t^  - 3t) /2 

p^(t)  - (35t^  - 30t^  + 3)/8 

Figui c 2 

Graphs  of  Legendre  Polynomials 


Our  Idea  Is  to  approxitiiatc  the  cc(|ucnce  of  each  coefficient  by  a func- 


tion of  time  whose  pararnoters  can  be  estimated  directly.  We  decided  to 


use  a function  of  Legendre  polynomials  for  this  purpose.  The  Legendre  poly- 
nofnlals  are  only  one  of  a wide  variety  of  functions  that  might  be  used  in  dif- 
ferent applications.  We  do  not  imply,  by  using  the  Legendre  polynomials  here, 
chat  they  are  necessarily  superior  to  Chebyshev  polynomials,  trigonometric 
polynomials  or  many  others  that  could  be  used  to  approximate  any  arbitrary 
function  of  time.  However  the  Legendre  polynomials  are  easy  to  visualize 
and  they  provide  an  adequate  expositional  device  to  illustrate  the  general  prin- 
cipals involved  In  our  technique. 

« 

Given  any  polynomial  of  time,  the  time  paths  of  the  coef f ic ients  o< j ^ and 


can  approximate  as 


li-c, 


( - 1 ^ t 1) 


fi.r  llo 


( - > ^ t ‘ » 


where  Pj(t)  is  the  polynomial  of  order  i at  time  C and  a^j  and  bjj  are  the  co- 
efficients of  the  ith-order  polynomial  for  the  intercept  and  slope  coefficients 
of  (1),  re-spoctlvuly.  In  this  case,  P is  the  set  of  I.egendre  polynomials  and 

the  first  five  of  these  are  given  in  Table  1 below  and  illustrated  in  Figure  1. 
IJnitE  of  time  are  chosen  so  that  -1  to  +1  spans  the  observed  number  of 
2 3/ 

natural  calendar  units- — . This  particular  structure  Is  used  because  the 

2a  / 

Legendre  polynomials  are  mutually  orthogonal  on  the  interval  [ -1,  1 
Thus,  each  polynomial  can  be  introduced  as  a separate  variable  in  the  es- 
timating regression  without  having  to  worry  about  mul ti-col li ncarl ty ; (which 
would  be  serious  problem  if  a non-ortliogonal  polynomial  vjcre  used). 

The  final  estimating  equation  is: 


+ b,,.R  + b,.tR  ^ + b,,  ((3t  - 1)/?]K  , -1 

OJ  ui,t  IJ  m,t  2j  m,t 
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For  n “ A,  the  coefficients  enn  evolve  in  ui*  to  a quaillt  fai.lilon  wl.li.h 
will  acl/iilt  a changing  trend  and  even  a limited  cyclic  pattern.  In  c.aih 
ease,  the  cr.tlnatecl  coe  f f lei  eiU  s of  the  h i (die  r- older  polyi.oml  nJs  n.lc.ht 
turn  out  to  be  Inslgnl f 1 can 1 1 y different  from  zero.  This  would  Imply 
Btntlonarl ty . 

Finally,  note  that  a test  for  the  existence  of  non-st.ition.nrl  ty  itself 
during  the  observed  sample  period  can  be  mode  very  easily.  The  comfioslte 
hypothesis,  = b^j  = ...  » 0,  using  all  except  the  zero-order  estimated 

coefficient  in  an  F-test , provides  a means  of  determining  tbe  probability 
that  1^  j In  (1)  (the  risk  coefficient)  has  actually  varied  in  some  way.  Si  mi  l.i  r 1 y , 
a j j = 32  j * • • . “ 0 is  a hypothesis  tb.at  neither  nor  h(  6 j^)  has  cdianged 

during  the  record. 


I 


.J 
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IV.  riff* Jf-1 '’"Ai. 


The  Data  Sample 


21' 


The  data  chnsen  for  1 11  ur.t  ra  t Ion  here  are  true  rnl<-,  of  letuin 

(Including  dividends,  if  any)  for  New  York  un  ' Americne,  SteeV  r.x- 

changc  listed  socutUie.-.  beginning  in  Inly,  l‘U,2.  Ob'e  i v,i  i i mis  were 

2_u/  , , 

available  daily  but  were  aggregated  to  weekly  Intor-al;  ' . Inc  last 
observations  were  for  July,  1969,  and  a total  of  365  wecKly  returns  were 

available  for  a large  number  of  r.c.  tiri  t i cs  . Since  there  were  jnute  limn 


adequate  data,  we  decided  to  use  only  a subset  of  securities  and  observa- 
tions, dividing  the  total  period  Into  subperiods  and  using  only  those  sU'CKv 
which  had  a full  record  of  obscrea f f ons . Furthermore , since  the  corput - t lons 
were  done  on  a commercial  computer,  wc  arbitrarily  decided  to  limit  the  major 
production  run  to  a fixed  nunber  of  minutes  of  CPU  time.  This  hal  the  effect 
of  eliminating  American  Stock  Fxch.ing.e  securities  v.'hosc  names  lanl.t'd  Ite  in 
alphabetical  order  since  these  were  positioned  l.ist  on  the  r.-apnciio  tape. 

The  final  sample  consists  of  930  securities  which  have  full  le  rids 
during  the  first  160  weeks  (from  July  5,  1962,  through  July  22,  I'  C'^).  Oi 
the  930  securities,  84  were  listed  on  the  Atiicrican  Stock  Kxchange  (aiiJ  the 

remainder  were  on  the  NYSh)  . There'  were  .some  stocks  e 1 im  i n.it  r ■ . tt'  1 u.d 

. 27/ 

spots  on  the  tape. — 


The  rcjason  for  choosing  160  weeks  was  that  the  robust  technique  requires 
the  sample  size  dlvideci  by  the  number  of  variables  to  be  an  integer.  Originally 
the  number  of  variables  was  2,4  and  10  so  that  the  sample  size  had  to  be  a mi:l- 
tiple  of  20.  Since  160  was  a reasonably  larger  mmibej  and  a rjultlplc  of  20,  and 
since  we  wanted  to  reserve  at  least  one-half  of  tlie  available  observal  ions  (,.io,' 


17 


365)  for  i>rccHctlcii  and  po'; t -i.ampU'  icstinj-,,  160  w.-i:;  dio-cn.  For  rt-asons 
to  be  discussed  shortly,  we  later  derided  to  do  the  robust  trlmiiin,  with 
the  two-variable  model  only,  so  we  could  have  chosen  a slifhtly  lary.er 
sample  (and  not  a multiple  of  4 and  10);  this  seerntd  hardly  woi  t h I hi> 
extra  programminp  effort,  so  we  did  r.oi  change  the  originally-chori  n sarplc 

size.  There  Is  no  reason  to  suspect  that  this  rather  hapliazard  technique 

? 8/ 

for  choosing  the  sample  had  any  Influence  on  the  results 

The  market  index  used  in  all  cases  v;as  the  Standard  6 Poor';.  nOO. 


An  Outli^  of  th ^_Re^u l_t s 

Since  the  results  comprise  a rather  large  amount  of  informatjen  on 
several  different  empirical  questions,  the  following  outline  will  serve 
as  a guide  to  the  most  interesting  part  for  each  reader. 


Section 
IV  A 


IV  B 
IV  C 
IV  D 
IV  E 

IV  F 


Containn 

Ordinary  least  squares  (01,S)  and  Fobust  (hOl.)  ve- 

resulls  for  individual  stocks  for  the  period  duly, 
1962  to  July,  1965. 

Tests  for  stationar itv  in  risk  coefficients  and 
other  parameters  for  the  sample  as  a vhole 
Analysis  of  individual  securities  \>ith  many  out- 
liers. 

Analysis  of  individual  securities  with  strongly 
varying  parameters. 

Tests  with  portfolios.  Results  for  the  first  160 
weeks  and  post-sam^plc  refitting. 

Tests  of  the  two-factor  n.odol  with  portfolios. 


fitting  may  give  better  results. 


we  will  piw.*.. .. 


A.  Ch.T  r (1C  t er  1 s r 1 c s uf  I iid  i v i ilun  1 Scon  j t 1 e s , July,  ! ■ !i  1 y , 1 S , 

Table  2 presentn  crosr.-f;cci  ioiial  KiatlsiicR  (aciect,  930  securities) 
for  regressions  calculated  with  tine  scries  for  each  security.  lor 
each  time  series  six  different  rec.risslons  v;ere  fit,  one  set  with 
ordinary  least  squares  COIS)  and  another  set  with  tlie  lolnrst  tech- 
nique, (KOB).  In  each  set,  tlie  tlirec  repressions  arc  'or  (1)  tiu* 
simple  Index  model,  (2)  tlic  model  with  linear  time  trends  in  the  co- 
efficients, and  (3)  the  model  iiith  up  to  a quarlic  polynomial  of  tine 
in  the  cocf  f icient.q ; (these  three  represent  model  (iDwith  n = 0,  1,  A). 


INDEX  MODELS  VITH  POSSIDLY  NON-STATIONARY  PARAMETERS  FOR  YORK  ANT) 
AMERICAN  EXCHANGE  SECL’RITIES^ , JULY,  1962-JL’LY,1965^ 


.963 


STATISTICS,  4-VARLABLE  MODEL 


'but  ^ome  mixlel  for  the  generating  process  was  needed  to  compute  large 


(IXTERCEPT  a,  OLS  -.200  .11 U -2.86 

QL’AXTIC)  R05  -.201  .822  -2.01 


The  tcfilc  pri'Tont  •;  tho  t.'il  cul  m cil  i rn  f f 1 c Icii  t ■ n ; l .iL  i 1 1 tr-  fii 

/s  * 

Hqj  and  from  r^juntion  11.  Ihoie  coc  ( f 1 c 1 prit  s r<i'i  (.  ■■(■nt:  tbc  pert- 

folio  risk  O^j)  .n-<’tcPpt.  J . the  .sii'pie  t t ^ "Ji' .lb ) 'V  .■ 

the  models  wltli  non- si  .it  i onnr  y ror  f f i c i on  t s , n ftud  b.  ^ are  the  cs- 
timated  values  of  the  pextfollo  lisk  and  the  Interccp'  oi!  average 
over  the  tiire  series. 


For  the  other  coefficients,  (a  andb  for  i>  0 ) onlv  thr 
t-statistics  arc  reported  in  order  to  save  space.  Since  tlio  entire 
interest  in  these  coc  f f i c i er.  t s it.  in  their  possible  deviations  f t oin 
zero  and  not  in  their  absolute  values,  icportinp  the  t-statistirs 
is  sufficient.  For  each  statist  ir,  the  cross-sectional  aiitln.'tit 
mean,  standard  deviation,  rinira^n,  fifth  perrentilc  , Tcdian,  93th 
percentile,  and  maxii.iurn  are  given. 

Serial  correlation  in  tlu'  residuals  was  Insignificant  for  all  irodols 
and  securities.  Over  the  .six  different  models  and  the  9-'i0  seonitiis,  tin? 
jninlniuni  computed  Durb  in'h'atson  statistic  was  1 .33  and  the  maxinn.'.  e.  is  ? . 7 , 
a range  that  can  bt-  attributed  to  chance.  The  me.an  ci  oss-sec  i i . na  1 vnlu.’s 
of  the  Durbin'h'atson  ranged  fto:,i  2.15  to  2.20  across  the  six  r.-dils.  T ni  ■ 
Indicates  a small  and  insignificant  negative  serial  dependence  in  the 
estimated  distu’'hauces , (with  .a  first-ordir  coriclntion  of  less  than  .1). 

TTie  major  differences  between  ordinary  least  sauarcs  and  robust  i eg- 
ression arc  evident  in  the  table.  As  mentioned  in  section  II,  the  robust 
technique  is  intended  to  guarantee  th.at  the  residuals,  ana  tiius  , he  i : ■ 
timated  coefficients,  follow  a Gaussian  prob.ahilitv  law.  The  success  C'f 
this  aim  is  measured  by  the  studentized  ranges  calculated  from  the  rcsi 
dais  of  each  model.  For  ordinary  least  squares,  the  entire  c ross-scc l ion.i  1 
distribution  of  the  studentized  range  lies  to  the  right  of  the  distribution 
calculated  from  the  robust  residuals. 

In  the  2-variable  model,  for  example,  tbc  95th  sample  percentile  stu- 
dcnllzcd  ranges  are  8.3''.  for  ordin.ary  least  squares  (01, S)  and  5.61  for 
robust  regression  (UOb)  respect  iv'^ly.  Tables  of  the  student  i.’cd  range 
show  that  the  95th  percentile  of  the  null  disttlbuticn  (Gaussian)  is 
near  5.75  for  a sample  size  of  160.  Tliis  is  quite  close  to  the  observed 
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Vfllue  for  tlic  roIai.M  cnsr  l>ut  Is  !.ir  l.olow  iho  obr.ctvcd  01, S vnlun. 

Kor  0I,S , cvfn  t!u'  nc-dlnn  obiifivad  r.t  ndcnt  i zed  i.mi’C'  is  far  at.ovi  I he 
93Ui  po  rcen  C i Ic  of  rlic  null  d i s t i I hu  i 1 <si  In  a ) 1 I In'?,  e r rp  i e.'  s i oi . J li  1 c 
Indicates,  of  course,  lli.it  the  obsei  ved  rcr.idu.ils  arc  quite  noii-C.iuss  i an 
when  ordinary  least  .''.quai  c;.  i .‘i  *n|ployed  vith  tliise  data.  Ihu'.  the 
conpnlcd  OI.S  t -i.  t a t i :.  t J c ;.  are  not  neressaiily  lellahle  mf.i'.ntes  of 
eigni f J cance , With  the  robust  technique,  however,  the  residu.ils  do  not 
violate  the  as.sumption  of  nornallty.  Thus  tests  of  s 1 r n i f i c.anc  e bared 
on  the  nortnallty  assurption  h.avc  face  validity.  Thi;;  is  very  ir.jiortanl 
In  testing  for  the  presence  of  non-statlonary  coefficients,  as  wesh.all 
sec  In  the  next  sub-section. 

One  other  fact  about  the  studcntiz.ed  range  should  be  noted:  Its, 

cr oss-sec 1 1 ona  1 standard  deviation  is  three  tin.-'s  larf.er  for  OI.S  than 
for  ROB.  This  is  no  doubt  due  to  some  extremely  large  values;  (the 
observed  maxitiiun  for  OLS  is  12.7),  which  indie. ite  gios.s  departures  from, 
normality  for  some  securities.  Whatever  the  source  of  these  departures, 
be  it  non-Gaussian  stable  probability  laws,  data  errors,  or  non - 
statlonaritics  in  the  parameters,  the  results  inaku  a strong  case  in  support 
of  the  use  of  the  robust  tcchnloue,  particularly  for  securities  with 
very  large  observed  values  of  the  sludentized  range,  (and  correspondingly 
large  potential  errors  in  estimated  parameters). 

The  explanatory  power  of  the  regression  seems  also  to  be  lirprove.i  bv 
the  robust  technique.  This  is  shown  by  the  observed  cross-sec L i or.al  dis- 
tiibution  of  R , which  is  shifted  to  the  right  from  the  OLS  distribution. 

Of  course,  this  is  hardly  surprising  .since  the  robust  technique  has  the 
effect  of  throwing  ovit  observations  that  do  not  conform  veil  to  the  es- 
timated regression  hvpcr-pl.ane.  One  should  be  vorv  careful  to  note  that 

this  supposedly  higher  explanatory  power  Is  really  just  a bias  in  the 

2 

R computed  by  the  robust  tcchnioue.  If  the  discarded  observation  had  been 

2 

^added  back  to  the  computation  of  R , (using  estimated  residuals  for  these 

observations  that  vjcre  calculated  with  the  robust  coefficients)  the  result 

would  be  a value  for  R^  much  closer  to  the  OLS  value.  For  the  same  reason,  t 

predictive  .ahilitv  of  the  ROB  ?quntions  will  not  be  better  than  the  OI.S 

2 

equations  to  the  extent  of  the  " Irprovcment"  In  R . This  is  because  the 
same  processes  which  generated  outlicr.s  during  the  fitted  sample  period  can 
generate  outliers  In  the  post-sample  period.  If  there  is  any  improvern  nt 


In  predictive'  nb  i 1 1 1 v I'/  poltij',  ( r<';^  (d.r.  ti'  I ' it , 
tbc  KOJl  ctx- f 1 1 c 1 u ' li.ive  Je'.;,  ( t ! r .i  I I eii  (iioi. 


; I V i i ] C'li  1 y be  t ( e .e  e;  I • 


Wc  b.cp.nn  tbc  !>tudy  by  ii'.ln,’  li.e  ifbe.t  netliod  te  i(  ;,eu  c.-etli  ■ ‘.oi 
each  model,  (2,  <'i  , and  10  vaii'ililes)  t;e]>at  aleJ  y.  lIi:  f o>  t ana : c 1 y , Ibia  is 
a very  t r.eac berou'.  oci-d'.u  e fi  i t!.f'  •!  .nui  10  le  r.odils,  t'S  •-.e 

soon  d i s cove  red.  Tliv’  rc  .a  son  s 1 1 : . t i * :■;  ' b o r t ht  y c>na  j p - 1 .■ik:"’  i - . 1 yr  ! 

as  rrgiessors  in  these  models.  Kec.ill  tli.it  the  robu'.t  t ethi;  i qiK  flr-t 
partitions  the  sample  Into  N/1.  snh-.saMpl  •.  s , nhere  b'  I',  thr?  ''.iC'ple  s i , 

and  k is  the  nurnher  of  vatlahlc  ..  A repression  i th.n  i . r , ’ ( -d  : ,<i  i ..s 
Bul)-samplc  and  the  nedi.ins  of  these  eoefflclents  at'’  t-ahen  as  tlu-  lir't- 
pass  eslinutes.  tin  f o r t una  t e 1 y , <.  .'l  a though  the  pal);.,  rial  • ere  at;;.  , .c.l 
over  the  full  time  span  of  160  ohsC' r va t i on.s  , tluv  ate  far  f ron.  c;  t'  ' ■ ’ ’ 

within  each  sub-  period. 

Tfiis  is  quite  appatC'iit  fi.aP.  a pi  line  ,u  tbe  p,i  . . !i  ' ti  e 1 i .p  : 

polynomials,  (Figure  1,  section  111).  ' f the  sue- r ..r-p  1 , . fot  tb.  :■  ,•  l 

method  arc  chosen  by  a natural  partition  along  the  Li'.r.c-  axis,  the  tntei- 

correlation  among,  the  pol  ynomi  a 1 ,s  within  eaeli  s ul,  - s ar  p 1 o '-.'ill  be  quite 

high.  For  the  10-vuriable  case,  the  axis  is  partitioned  into  ICO.ijO-  16 
equal  sub-sec  t i on.s . A one-slxtecnt  b sub-section  of  tbe  time  axis,  rhot.e;'.  an 
anyvb.cre  at  all,  v.’i  1 1 result  in  extte;  ly  In' pb  i nt  c r- depend.:  ne  i es  .nv. 
the  polynomials.  The  6-varl.ible  ca;.e  i.s  not  ns  tt  .uhli.ig  1 c:  .;u  .c  . 

first-order  polynomial  and  the  first-  order  polynoni.al  nulliplie.l  ly  the 
market  return  are  used.  ilowever  , the'  tuib- s.e  c t i ons  arc  sm.aljC!  , r.ivh 
being  1/AOth  of  the  total  sample.  Thus,  some  neatly  sinp.ular  matrices 
have  a high  probability  of  occurenre,  ;. imply  becau,..'  Icuir  suceessi-.'e 
market  returns  of  approxii.mt  ely  the  tune  r.aguitn.ie  .u  e i.ithet  1'.  ■ -V  to 
occur  at  lca.=  ' once. 

We  tried  to'  remedy  lliis  induced  uiul  t i-col  1 incari  t y first  by  choos  1 ng 
sub-samples  by  some  ' ’thod  othci  than  the  natural  ordering  along  the  time 
axis.  For  example,  we  tried  using  observations  numbered  1,  11,  21,  ....  151 
for  the  first  sub-sample,  2,  12,  22,  ...  for  tbe  second  sub-sampic  etc..  In 
the  10-variable  case  and  1,  5,  9,  etc.  for  tbc  ^.-variable  regressions.  This 
alleviated  considerably  the  mu  1 1 i -co 1 1 inca r i t v problem  for  the  lO-variahle 
case  but  it  did  not  help  tlie  ^-vari.ible  case  at  all,  and  there  ^.’cre  still 
some  sub-samplif!  in  botli  cases  wbosi;  moment  matrices  had  cxtrcn.ely  low 
determinants. 


Finally,  vc  to  do  the  first  part  of  it.i-  lolmsl  tc(  1ml  C(i'r  , til.’ 

detormlnatlon  of  outlier  obr.oi  va  1 1 ons  , with  the  tvo- var  iabic  rnodcl  alcuic. 
Then  tlic  second  pa-;-;,  vhi  cli  consists  of  cohtj-ui t a t i in-  ordinary  least  sciti.irts 
on  non-excludcd  observations,  was  done  for  the  2,  I*,  aiul  10 farialle  model 
but  using  the  observations  determined  hy  the  ?.-vaiiablc  first  pass.  Ad- 
mittedly, this  is  a somewhat  atblliaiv  procedure  an!  o.m  would  exptet  a 
priori  that  it  would  undorslatc  the  true  worth  of  the  robust  iechnie,uG 
relative  to  OLS  for  the  A-  and  10-variablc  models.  However,  the  special 
nature  of  the  independent  variables  In  these  models  made  it  necessary 
and  we  can  do  no  better  than  to  keep  it  in  wind  wl.en  interjir  ct  ing  the  i e- 
sults . 


i 
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B,  Test.liir,  f tl.«~  I'  <'.  s l V nc  c ti' 'i-.- s t J t 1 on3  r y ;»  i a "i-l  < !'■; : A IMrS'  !'J 

One  of  ihc  basic  «l  r rac  t i ■ iis  of  on  nctlio'i  Is  its  c.i,  ability  t('  provlHo 
a statistical  test  for  rhf  e-x  la  l one  e <1  i.  .n- s t a t i onn  t i i / . Uc  1 1 1 ns  I i a;  f 
this  by  tests  of  stationarlty  of  the  c f f j c 1 one  s of  (1).  In  pa  r 1 1 c.i  1 o r , 
the  stability  of  risk  parai-.cter  ( p j in  (1))  over  tj.p  i ■,  a n:,'i!tcr  of 
great  Importance  to  portfolio  managers  arui  other  in-’estors.  They  would 
like  to  be  able  to  use  h 1 s t oi  i c a 1 1 y- es  i i t a t ed  risk  par.-ntter!.  for  cuirei.t 
portfolio  selection  and  would  also  like  to  be  able  to  extrapolate  li.t' 
safely.  In  addition,  the  intercept  leincXj  is  somet  i'-er  i nr  er]- r e t eri  as 
the  "extraordinary"  return  to  a stock,  above  and  beye  ud  the  ;i  rr.ia  1 co:;,pen- 
sation  for  its  risk.  Thus,  it  too  is  imporlnnt  in  so.ie  contexts.  We  will 

report  evidence  on  ti'.e  .>tat  It-nari  ty  of  a l)oth  j-a  ran  c t e i s . 

The  central  problem  in  testing  for  tlie  over-all  existence  of  non- s La  t ! r ,i  i r 1 r y 
1b  the  a priori  lark  of  Vtncnvledgc  of  the  direction  of  rh.inge.  Dut  ing  the 
sample  period,  we  might  be  lucky  ami  obr.eivc  rlsl  coefficients  on  all 
securities  drifting  together.  Certain  reasons  foi  dtift,  siicli  as  c lianc.i  r.;; 
BltitudcG  toward  risk  on  the  part  of  all  inver.to’.s,  airi.L  leat.onil  !y  1 r 
supposed  to  cause  such  a gcner.il  and  concurrent  mo.  emen!  . H this  li ) .d  e.rv.;-, 
ve  would  observe  that  ‘.oinc  or  all  of  the  coefficients  a-'-sociated  with  jcly- 
nomials  of  time  in  the  cstin'.ai  ing  equation  (|  | ) were  .signi  f icant  ly  r.  - .(to 
on  average. 

Other  reasons  for  uon-s  t a I i onar  i t y , such  as  cltangc".  in  individual  f i t c . r.ij  ii.al 
Etructures,  would  produce  non- cone  ui  ten  t movonentf.  in  risk  coefficient'..  Pir. 
coefficients  might  increase  and  othei.s  decrease,  depcitding  on  the  specittc 
circumstances  of  each  firm.  Since  the  average  value  of  the  lisk  cot f f i c i-o.t 
must  be  very  close  to  unity,  these  non- concur  ren  t non- st  .i  t i on.ii  i t i e s would 
not  cause  the  mean  cross-sectional  values  of  the  polynomial  roofficicais  of 
(11)  to  deviate  from  zero.  They  v/onld,  however,  cause  the  estimate!  polvnonla! 
coefficients  to  follow  a different  c ro;;s-soc  l i ona  1 d f s 1 1 i t-u  t itnt  tlian  th.it 
expected  under  the  null  hypothesis  of  non-s  I a I i ona  t i t y . The  distiibution  mip.ht 
very  well  be  located  at  zeto  hut  the  (ro.sr.-section.il  diipeision  would  be 
larger  than  the  anticipated  null  dlr.peision. 

This  Implies  tli.it  tests  for  non-concurrent  non-r.  t at  ion.ii  i t y are  highly 
dependent  on  a knowledge  of  the  distribution  of  the  estimates  from  (II). 


contrast,  the  robust  rej’ res;,  i on  t<'c!.niqiJo  f n.n  anr  t-os  asyi'qit  tit  Ic  nor:  .a  1 i ty  . 
Furthermore,  studenti/od  ranges  repotted  In  Table  2 support  the  tonlii't 
that  our  chosen  sample  size  Is  sufficient  ^o  rale  valid  an  asrur.pt  . on 
of  Gaussian  disturbances  from  the  rcdiust  ii.odels.  It-eMfoii  , oc  will  n ' 
the  robust  results  in  testing  for  nOn-stationat 1 ty  since  we  ran  be  i tasotial  ly 
confident  about  the  sampling  distribution  of  the  t-iat!os  ftt  this  r e rhi.i  y . 

Under  the  null  hypothesis  of  totally  stationary  ccef f Icients , the  higher- 
order  rohus  t -cs  t Ir.ia  t ed  toe  f f 1 c 1 cat  s (11),  ^Ij’ 

normally  distributed  with  itiean  zero  .md  i t .ind  1 1 d di-v  i a i 1 1 ir,  c. ' 1 1 i.i  t ■ : I y tin  : i 
computed  standard  errors.  Tlie  t-ratios  are  disltlbutcd  acroiJlng  to  the 
Student  lav/  with  about  146  and  140  degi  ces  of  freedom  for  the  4 and  10  variable 
models  respectively.  (I'cgrecs  of  freedom  differ  across  securities  according 
to  the  number  of  rejected  outliers)  . This  Student  ^ ^v;  Is  verv  close,  but 
not  exactly  equal,  to  the  st  nndard.irized  Gaussian  l.iw.  For  d.f.  ^ 120,  the 
90  pur  cent  interf  racti  le  rangt'  is  3.316  for  Student  '.'hile  it  is  d.290  1.; 
Etandardarl  zed  Gaussian.  The  difference  is  snaller  for  all  but  one  of  the  0'’^' 
securities  since  the  minimum  sample  size  is  128  and  the  next  lowest  .‘wir-plc  size 
Is  134.  Therefore,  using  the  Gaussian  approximation  will  result  in  only  a 
trivial  numerical  bias. 

Tlio  90  percent  range  is  itself  asymptotically  normal  with  standard  dcviat.'on 
approxlm.a tely  equal  to 

Oj^  9-(l/f  p^)(2(.045)/N)^^^  ■=  0.234 

where  f.05  is  the  standard  Gaussian  ordinate  at  the  5lh  pciccntilc  and  N is 

29/ 

the  sample  size  (which  we  have  taken  equal  to  155  in  this  calculation) — 

Confideitce  regions  for  t^c  ranges  of  estir\aied  coefficients  can  he  construettd 
with  these  numbers.  For  example,  tlic  range,  3.200  + ( . 2 34 ) ( I . 96)  , oi  [ 2.83,  l.t' 
is  the  95  per  cent  acceptance  region  of  the  observed  90  percentile  range  for  the 
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null  hypothesis  of  s t .1 1 i oti.ir  i t y . It  tlu‘  ohsc  I'viil  v.iiif.i-  foil;,  within  fliis 
interval,  the  null  hypothesis  cannot  he  rejected.  Table  ■ presents  test 
Statistics  associated  witli  this  idea.  The  table  gives  computed  values  of  a 
"standardized"  90  percentile  range,  z,  defined  as 

t s (R.9  - 3. 290)/. 234 

where  R.9  Is  the  ef.tlir.aied  90  percentile  range  (computed  a;;  the  difference 
between  the  93th  and  5th  cross-sectional  percentiles  of  the  t-ratios  of 
table  2).  Values  of  z arc  starred  if  they  imply  an  observed  range  outside 
the  95  percent  acceptance  interval  and  thereby  reject  the  hypothesis  of 
stationarity . 

Tlic  table  also  presents  tests  labelled  t/  (o^/>^)  for  the  significance  of 
the  cross-sectional  uiean  t-ratio  and  c/o^/A N)  for  the  significance  of  the 
cross-sectional  mean  coefficient.  The  symbol  t denotes  the  arithmetic  mean 
t-ratio  in  Tabic  2 and  0^  denotes  the  ci oss-scct ion  standard  deviation  of 
the  observed  t-ratios  . c is  the  arithmetic  cross-sectional  moan  of  the 
estimated  coefficient  and  0^  is  the  cross-sectional  slnndnid  deviation. 
Values  significantly  different  from  zero  at  the  95Z  level  are  starred. 

In  comparing  these  tests  with  the  range  test,  one  can  think  of  the  1 .ingc 
as  measuring  non-concurrent  changes  (across  securities)  in  the  model  s 
parameters  while  the  mean  t-ratio  and  mean  coefficient  measure  concurrent 
changes. 


In  both  the  4-variable  and  10-variable  models,  there  seems  to  be  evidence 
of  significant  non-stationarity  of  both  types.  For  example,  the  "slope 
trend"  (associated  with  the  first-order  polynomial  coefficient  in 

(11)),  has  a t-ratio  which  is,  on  average,  significantly  negative.  The 
10-variable  results  indicate  that  the  averj^^drift  in  the  risk  coefficient 

was  not  a simple  linear  function  of  time  (as  it  must  be  regarded  by  the  4 
variable  model  which  can  accomodate  only  the  first-order  polynomial).  Both 
the  quadratic  and  the  quartic  polynomials  .rssociated  with  slope  changes  had 
significant  mean  t-valucs  in  the  10-v.iri.able  model.  In  fact,  the  cross- 
sectional  mean  slope  trend  is  only  marginally  significant  in  the  10-variablc 


TARI E 3 


Tests  for  the  Presence  of  Non-Stat  ion.iry  I’aranictcrs  Ajiiong  930  NVSI'. 
and  AMEX  Securities,  1962-65 

(Robust  Kegi  c ;;s  ion) 


Est  itnated 
Parameter 

No.  of 
VKl.S.  in 
model 

Standa rdard i zed  .9 
range  (z) 

t/((J^/,^J) 

c / ( 0 / »^' ) 
c 

Intercept 

Trend, 

4 

-.769 

7.90* 

9.09* 

Slope 
Trend,  6^ 

4 

3.54* 

-3.56* 

-2.65* 

★ 

Intercept 
Trend,  2^ 

10 

-1.07 

8.93 

9.74 

Intercept 
Quadratic,  32 

10 

-1.41 

.725 

1.91 

Intercept 
Cubic,  a^ 

10 

-2.35* 

-6.29* 

-5.63* 

Intercept 

Quartic, 

10 

-2.35* 

-7.45* 

-5.69* 

Slope 
Trend,  6^ 

10 

1 .83 

-2.22* 

* 

-1  .^.7 

* 

Slope 

Quadratic,  62 

10 

1.15 

9.18 

7.94 

Slope 
Cubic,  6^ 

10 

.0427 

1 .80 

1.49 

Slope 

Quartic, 

10 

.940 

-4.22* 

-3.19* 

^Rejection  of  the  hypothesis  of  stalionarity  <it  tlic  95  per  cent 
level  of  significance. 
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model,  wliich  means  iliat  most  of  its  variation  was  iion-l  i near , 

The  mean  intercept  also  sccias  to  have  changed  s I gn i f leant  1 y (and  non- 
lincarly)  during  tin-  sample  period.  This;  is  of  com  sc  not  too  sur[>r  i r.i  up. 
since  the  intercept  should  he  a n<f..it  ivc  linear  function  of  the  slope  acc- 
ording to  some  theories.  Indeed,  the  first-order  polynoniial  coefficients 
arc  in  agreement,  since  they  arc  positive  for  tlic  intercept  change  and  ne- 
gative for  the  slope  change.  They  arc  more  significant  for  the  intercept 
than  for  the  slope,  however,  and  this  could  indicate  a non-stat ionary  risk- 
less rate  of  interest  or  the  niovcrncnt  of  some  other  omitted  factor. 

The  time  paths  implied  by  these  mean  coefficients  are  given  in  Figure  4. 

This  plot  simply  assumes  that  mean  estimated  coefficients  of  the  higher- 
order  polynomials  accurately  depict  the  actual  paths  of  the  slope  and  iiUui- 
cept  during  the  sample  period.  We  must  emphasize  th.ii  this  graph  is  biased 
toward  showing  no  movement.  If  Che  index  used  as  the  explanatory  v.ni.Tulc 

had  been  composed  of  the  storks  in  the  sample,  the  mean  eross-seet  ion.al  f. . 

J • '■ 

would  have  been  unity  for  every  period,  by  const  rue  l ion . Ttius  it  could  iiot 
have  been  non-stationary . Figure  2 shows  movement  only  because  the  index  was 
not  composed  of  the  sample  of  securities  weighted  in  the  same  proportion. 

The  tests  for  non-concurrent  changes  also  find  a few  significant  values. 

The  slope  trend  in  the  4-variable  model,  for  example,  has  a larger  90  percent 
range  than  one  should  have  expected  under  the  hypotlicsis  that  all  coefficients 
are  stationary.  This  agrees  witli  the  first  two  polynomials  froni  the  10-variable 
model.  (They  are  significant  together  but  not  individually).  Evidently,  in 
addition  to  many  securities  displaying  concurrent  drift,  there  were  also  some 
which  drifted  significantly  in  a direction  opposite  to  the  aver, ago.  This 
implies  that  the  concurrent  drift  was  actually  more  substantial  Lliaii  the 
average  coefficient  indicates  since  the  average  coefficient  contains  some 
securities  which  moved  in  opposition. 

There  is  unfortun.itcly  a statistical  problem  with  these  tests:  they 

assume  that  the  cross-sect ion.il  distribution  of  estimated  coefficients  con- 
stitutes a random  sample.  If  there  is  i nt ordepcndouce  among  the  estimates, 
both  tests  will  be  biased  (but  in  different  directions).  For  example,  if  there  i' 
positive  dependence  among  the  cstim.itcs,  the  computed  stand. ird  error  of  the  nxau  i 
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v<iluc  (6|^//N)  will  I'O  lo'wcr  tliat^  lln  f.i  .itul.ii  .1  ciioi  lti.it  ~.\ild  have  liccii 
obsctVL'cl  ill  fl  j;criiiiui'  i .milom  saiin'lc.  Ttiir.  v.'i  1 1 ovci-l.il''  i lif  i i',n  i f u i nr  i‘ 
of  tlic  mean  t-valuo  (ur  .all  r.cc  in  i I i (.nnJ  llm;.  will  i lir  ((immiiut 

drift  ajipcar  r:ore  s i p.n  i f i r ant  tl..in  it  trally  is).  On  l lie  otlui  li.ind  , tlie 
Staiidard  error  of  the  t.anplc  .9  raup.o  does  not  depend  on  the  observe  t i rr.-.r. . 
(It  is  computed  from  the  known  distribution  of  a random  sample  of  ordci 
statistics  from  a Gaussian  law).  Ihorefore,  positive  deptndcnce  ais'  np, 
the  estimates  would  cause  this  test  statistic  to  under st  ate  the  true 
significance  of  non-conru r ren t movements  (because  the  ohseived  range 
would  be  smaller  than  a tango  obtained  from  a random  sample).  The  actual 
interdependence,  if  any,  would  depend  on  the  cross-sectional  dependence  in 
the  disturbances.  King  tl96G  J found  a very  small  degree  of  pot.itive  crOsS- 
section  dependence  in  the  di  sturbancis  from  a versicn  ot  the  ? v.iria.hle 
model  but  even  a small  degree  of  covaiiation  can  cause  a latgi-  bias  in  tbs'* 
estimated  significance  of  the  mean  coc f f i cicnen t s or  the  mean  t-ratios. 
However,  since  the  prc.scncc  of  interdriiendcnce  would  bia'-.  the  raiipp’  ii  i 
in  an  opposite  direction,  and  since  both  types  of  tests  indicate  sos.e  Sig- 
nificant non-stationarit  i os  in  the.  coefficients,  we  arc  safe  in  concluding 
that  some  kind  of  non-st a t i onar i t y was  present  duritig  the  sample  per.nd, 
even  though  ve  should  he  hesit.int  to  st.ite  that  it  w.is  concurrent  rather 
than  non-concurrent.  In  Section  IV-K,  this  question  will  ht’  re-cxe.nii  lu  d 
with  estimates  obtained  fro.si  portfolio  returns,  which  \.i  1 1 be  used  to 
obtain  an  estimate  of  the  cross-sectional  interdependence. 
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C , PcfcctinK  PniiS'.ial  Ohsorvfi  r i ons  | 

Another  potentially  luiefol  capability  of  the  ttclmlqvif  It.  the 
Identification  of  unusual  cases.  UTitn  the  same  model  is  applied  to  a 
large  number  of  cases,  an  autcanatlc  method  for  detecting  departures 
from  the  norm  would  be  highly  vnlual>le.  F('r  eKample,  a l ank  r.ii  t,ht 
depict  the  credit  risk,  of  a large  number  of  its  term  loan  (usLoner;.  by 
a model  using  the  periodic  accouiUing  information  of  each  client.  Sud- 
den and  unusual  changes  in  the  iv.odel's  fit  could  indicate  changes  in 
credit  worthiness. 

Our  technique  provides  two  Indications  of  such  "umisual"  circum- 
stances. The  first,  to  be  derived  in  this  section,  relies  on  observa- 
tions detected  by  the  robust  procedure  as  lying  outside  the  non.ial 
range.  The  second,  to  be  discussed  in  the  subsequent  section,  involve.s 

unusual  movements  in  the  estimated  cocfflcit-nts  of  the  linear  mi'del.  - 

Again  reverting  to  our  example  model  (I)  with  slt'ck  return  data, 
the  cross-sect  iona  1 d i .s t ribut  1 on  of  the  number  of  uon-exc  ludt-d  obser- 
vations from  the  first-pass  robust  regression  is  given  in  Table  tt . 

Table  ^ 

The  number  of  non-oxc  1 uded  ohsorv.n  t ions  ^ us  i ng  tbc  I 

robust  method  on  930  NV.Ib  and  A.'IFX  .sccur  i t i c;. , >j 

July,  1963  - July,  1965  '1 


• 

Mean 

Standard 
devia  t i on 

Mlnimim 

5th 

percent  lie 

median 

95  th 

percent i 1 o 

^M.TxirAT 

j 

• ■ 

155. 

3.50 

128. 

149. 

156. 

160. 

1 t-O  ! 

1 

Maximum  possible  is  160. 
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Klpurcs  3 and  6 provide  furtbor  information  alout  tbc  effect  of  robust 
regression.  Figure  3 gives  the  cro;.s-sec t I onal  frequency  d 1 s I i i bu t inn  of  tbc 
ntimhor  of  oiitllrrs  (which  Is  160  ininus  the  nunbet  of  included  obser  v,i  i ions). 


930  NYSE  and  AMEX  nccurlties,  1962-65 
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This  shov.’s  that  58  of  the  930  securities  had  no  excluded  ohscr  va  t Ions . 


For  these  hccu r 1 1 it- s , o(  rourne,  l h(  1)1, h uiiil  lolmst  i".I  liinites  v/ei 
Identical.  The  moat  frequent  number  of  outllero  wa;.  In  tlie  r.inf'.e  2 to 
4,  which  is  a very  small  percentage  of  the  raniple.  There  were  some 
eecuritles,  however,  with  a significant  proportion  of  excluded  observations. 
These  will  be  the  subject  of  further  analysis  in  sub-section  C. 

Figure  ^ presents  the  number  of  outliers  as  a function  of  time.  The 
Interest  in  this  chart  derives  from  jiossible  non-s ta t iona r i t 1 cs  in  the 
dispersions  of  Individual  security  disturbance  terms.  If  such  non- 
statlonar it ies  are  interrelated  across  sccutlties,  v;e  should  observe  a 
pattern  in  the  time  path  of  fhe  nun, her  of  outliers.  Indeed,  Figure^ 
shows  that  the  five  weeks  v;ith  the  most  outliers  occured  before  v.’eek  30 
(which  ended  January  24,  1963).  There  is  also  some  serial  dependence  ever 
time  in  the  number  of  outliers.  This  would  be  consistent  with,  for  example, 
non-stationary  but  slowly  varying  dispersions  of  disturbances  for  all 
securities. 


To  check  whether  part  of  this  is  due  to  unusual  moverionts  In  the  S&F  500 

Index,  caused  by  abnormal  changes  In  one  or  more  of  Its  component  securities, 

30  / 

we  calculated  the  follovjing  relation  between  the  natural  logarithm — of  the 

number  of  outliers,  N,  and  the  current  and  lagged  absolute  values  of  the 

Index  return  R , (The  absolute  value  of  the  market  index  return  is  a proxy 
' m ' 

for  Its  dispersion) : 


Log  (N J = 3.05  + 5.81  )r  I + 9.14  Ir  I 
® ^ (50.1)  (1.39)*  ' (2.20) 


t=2,...160. 


Durbin-Watson  1.15;  Studentized  Range  5.79 
The  numbers  in  parentheses  are  t-  ratios. 

Because  of  the  low  Durbin-Watson,  the  regression  is  mls-spec 1 f led . As  a 


i 
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rough  remedy,  we  added  a lagged  dopendent  variable  as  a regressor. 


This  resulted  i ri 


log^(N^) 


(7.82) 


d 3.96 
(1.05) 


+ 5.66 
(1.A9) 


K 


n , t*- 1 


+ .625 
(6.08) 


log 


t-2.. 


.160 


DW  “ 2.13  , SR  » 5.79 


The  significance  of  the  lagged  value  of  the  nur.bor  of  outliers  justifies 
our  previous  staternent  about  serial  dependence.  This  regression  shov.’s  also 
that  the  influence  of  the  market  on  the  number  of  outliers  Is  not  very  strong, 
which  implies  In  turn  that  most  of  tlic  outliers  arc  due  to  events  specific 
to  individual  securities  or  specific  to  groups  of  sccuritli’s  In  the  sample  and 
are  not  due  to  unusual  changes  in  the  general  market  index.  Note  well  that  tb. 
average  value  of  R is  ahovit  .008  vUile  the  averni'C  v.nluo  of  lug  (N.)  is  about 

in , t c t 

1.6.  Thus,  even  though  the  coefficients  of  1r  I and  IR  ,1  are  larger  than  t!ie 

'm,t'  'm,t-l' 

coefficient  of  the  lagged  value  of  the  number  of  outliers,  their  average 
affect  on  the  dependent  variable  is  only  about  l/20th  as  great. 


T.io.H,  v.f  havo  Jrtiiimsi  rai  fil  UiM  ih'.*  lohu-t  H'ClmKjt.r  f' i i ••  i no  i c-  , 


r 

! 


! 


r 
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by  H I /I  I io  , r - pi*t  if  1 od  r"- 1 I «'  r>li  s ••r-  ..  i i - - vl.  1 « I*  *<  r c 

not  conf.ruoni  with  the  majority  t’f  otj-.c  i va  i l ons  . lii'*  acii.al  nuinl.rr 
of  "outlier'"  varle*.  a^.r«>'.s  n<’CuMl:‘  Vo  a «yt"f.l, 

part  of  tilt'  ^ a r I a t i ' III  n-M  bt-  ,i*.crlb''i  i ■>  rbaoce,  f-ln'o 

ra/ifliw.i  iM  / f I'j  rnro!i  .ii.iOi.f;  j li/>  9J0  ;.>  < m i i : nr  r bi.i;nd  !■<  tn  r.  ; . Jt  j-,  .il*i 

I r ue  , I , i ). .1 1 i.»  t II I 1 r 1 « ; i.  i t l.  .i  1 .i  i j,o  p i - • pi'i  i i • u ■ I i . i c > i i i • ! ■ • i .m 

I I « 'li«i  M 1 t i<  • > 1 t 1 i I • I / I ' 1 1|  •>!  b<  1 ■ > . i>  1 1 1 I ( I ••  ii.i  v«  < . I 1 1 1 < I.  • Ml 

am!  in  I ( (|  lit  <'  v<  . M.  w 1 I b p i r > t • i t t > . u > ■'  i { 1 1>>‘  ( I < , I • | • i i ■ I , ' - i .mi 

effort  to  I'l'tiiliy  sui-h  uniqup  c i rt  u ■ flue  r , wo  bjve  • ni  at  ••'i 
Out  for  detailed  analysis  all  seriir  i t i »•>  with  irore  than  JO  out- 
liers (O'jt  yf  l,,i^  t>bsf  rva  r 1 oils)  durinj^  the  sai.pli'  period.  Tiore 
were  ^*8  socurltics  in  this  class,  or  slib^i^y  r' i e than  five  pt-rcont 
of  all  securities.  They  are  p.iveu  by  nane  in  Table  5,  along  with  the 
nurbors  of  oitlleri  t>Sr-erv'cd  for  each  (In  any  other  application, 

of  course,  the  peicfiiinge  of  cases  to  he  singled  out  for  f it 'her  ana- 
lysis could  hu  diffen-nt.  Here,  we  arc  merely  1 1 1 ;is  t r a c i n,.,  fhe  seni- 
autonatic  nature  of  the  culling  ptoccss). 

Tlie  suh-sa-nj;  Ic  consisting  of  tlu’se  48  h i gb.- ■ 1 1 1 i • r sf>ick'  diff't- 
in  some  Important  respects  from  the  total  sai',plc.  Son-  »!  i j f ei  ein-es  coilJ 


f ' 
1 ■ 


have  been  foreseeable.  For  example,  the  mean  siu«>eiitl:M  ' range  of  le'itJuals 
cf»nputed  by  nr<llnary  le.ist  t,q  iarcs  (01. are  7.22,  hiu'  . .21  foi  tin  / , 

4 and  10  viirlnble  r.vidfls  respectively.  This  is  O.ft?  I»i|!fi  t'^iu  fh<'  'll.  r.’  a i 
Bt  iiden  I i /C’H  langcs  for  all  securities,  (See  I.tMc  2),  v!ii.h  itplii-. 
disturbancis  fc’r  ilie.se  h. ) git -on  t 1 i o r securities  dcvi.ito  fuithei  lion  (..u,  r.  i . 
than  do  all  securities.  After  application  of  the  tcdu'.t  mcMi'  i,  the  r '• 

Ftudcnti7cd  r.ingcs  acro>-.s  the  4ft  .stmiities  drc'p  to  b.Cft,  !<.08  nnd  S.16. 

These  are  vi*ry  close  and  even  slijlitlv  less  than  the  rii  an  s t uJcii t i r uH  i.-inys 
computed  with  the  rolmst  methed  on  all  930  sccuiitie'-.  Thus,  even  thcup.h  the 
disturbances  seen  to  h.we  been  very  non- Gaussian  over  tie  full  sample  record 
and  more  noj»- Cau.s.s  inn  for  the  higlr  ciitllcr  group  than  all  .sncuTitirs,  thi 

irjectlon  of  a re  1 .i  t i vo  1 j'  laipe  nu'’!r’r  of  oh  si- r va  l i <-ns  ‘.eims  to  h.'.ve  tuoj,!it 
the  derlrc'd  irsult  of  Cju..l.m  d i r- 1 r i bu  t cd  cMim.ilc;. 

The  hif.h-outlier  sul--s.n:iph  jl;,,,  uiflc,-  fiom  tie  snplc  of  .ill  scf.it  jLJ.^ 

Ill  ways  that  could  not  have  hei  n f o i cset’O . The  aviiape  prcpoitjon  of 
cxpl  0 i ncd.  v.ir  J a t i on  , (H  ),  is  nuch  Joi/ri  foi  these  .securities,  hriiip  lose 
than  .07  for  01, 5 and  about  .08  for  This  Is  only  balf  the  prepoitlcn 
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TABLV;  5 

Securities  with  >'..riiy  UnuMi.il  Obsf  i vn  t 1 ons  , July,  lOlwl-Julv,  19Gb 

Kfw  York  am!  Anuilcai\  Kxcliaoj’.i'a 


NAME 

NO.  OF 

ouiiMi:!-:: 

l! 

1(  NAME 

i ' 

NO.  OF 
OHTLl! KS 

Acme  Precision  Products* 

1? 

1 

1 CfuciAl  I’orcland  Ctment 

13 

Allied  Control  Co.* 

13 

1 Ginn  A Co. 

lb 

Araerace  Corp. 

12 

Gordon  Jewelry* 

14 

American  Bakeries 

11 

Great  Basins  Petroleum* 

17 

Ancorp  Nat'l  Services 

13 

Cr(Mt  Nortlieia  lion  Ore 

15 

Argus,  Inc.* 

18 

! Ilartfied  Eodys,  Inc.* 

11 

Atlas  Corp 

19 

Indiana  General  Corp. 

14 

Baltimore  & Ohio  RK 

lA 

Libby,  McNeil  A Libby 

16 

Bausch  & Lomb 

11 

"crit  t-Chapman 

32 

Beatrice  Poods 

11 

Minn  Lntcrprl.ses 

14 

i 

Bell  Intercontinental 

lb 

Monon  RR  (B) 

lb 

1 
, f 

Briggs  Mfg. 

13 

J.J.  Ncvd'crry 

11 

City  Stores 

17 

Oklahoma  Cas  and  Elec 

13 

t ‘ 

Colt  Industries 

14 

Pacific  Tin 

lb 

Continental  Materials* 

11 

Random  House,  Inc. 

11 

Countryv.'idc  Realty* 

lb 

Rapid  American 

18 

1 

' t 

Dorr-Oliver,  Inc. 

11 

Reliable  Stores 

lb 

Electronic  Commun.* 

11 

Sland/ird  Oil  of  Ohio 

1 3 

4 

El  Paso  Natural  Gas 

11 

Stone  A Webster 

12 

Esquire , Inc . 

11 

Texas  Pacif i c Land 

11 

Family  Finance 

11 

Tractor  Supply  Co. 

12 

i 

Federal -Mogul 

11 

Tianswcstein  Pipe  Line 

13 

Fischbach  & Moore 

12 

I'dylltc  Corp 

12 

Forest  City  Fntrepriscs* 

11 

We  lb  i 1 1 Cor p . 

26 

J 


*Anicrican  Kxclianpe  Listed  Security 
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of  explained  veriation  for  all  fici.  in  i t i rs  . 

Tabic  6 presents  evidence  on  still  anoctier  significant  difference,  tlic 
absolute  magnitude  of  the  mean  portfolio  risk  coefficient.  for  this 
48-sccurity  s ub-s ample  , the  mean  risk  coefficients  ai o s i gn i f i can ! 1 y lower 
than  for  all  securities  (cf.  Table  2).  I’urtlierinorc , the  estimates  obtained 
by  the  robust  method  are  si gni f i cant  1 v less  in  absolute  magnitude  than  those 
obtained  from  OhS . This  is  not  matchod  incidentally,  by  lower  t-values 
for  the  robust  coefficients--^.  In  fact,  the  robust  estimates  are  slightly 
more  significantly  different  from  zero  than  the  OhS  estimates  although 
both  arc  less  significant  than  the  average  coefficients  for  all  securities. 

TABl.K  6 

Estimated  Bisk  Coefficients'*'  for  the  48  High-Outlier  Securities 


NYSE  aud  /d-IKX,  1962-1965 


Model,  (No. 
of  variables) 

Regression 

Teciiniquc 

Mean 

Standard 
Devi  a t ion 

Minir.um 

5th  I'cr- 
ceiu  i 1 e 

Med  i .an 

95th  Per- 
cen  tile 

M^x  irui:" 

_ 

2 

OhS 

.82,1 

.438 

-1  .04 

.0918 

.904 

1.33 

1.64 

ROB 

.596 

.345 

-.502 

-.219 

.627 

] .00 

1.43 

4 

OLS 

.808 

.451 

-1.08 

.0598 

.819 

1.33 

1 . 58 

ROB 

.603 

.344 

-.501 

-.0369 

.606 

1.05 

1.35 

10 

OhS 

.743 

.540 

-1.86 

.0478 

.772 

1.45 

1.54 

ROB 

.564 

. 365 

- .929 

■ .00750 

.573 

1.00 

1 . 36 

* l.e.,  Estimated  cocf f icienls,  h . iBc  zero-order  polynomial  in 

model  (11) 

1 

1 

1 

1 

This  is  somewhat  puzzling  result.  for  a group  of  securities  with  relatively 
large  numbers  of  unusual  observations,  we  miglit  have  guessed  that  KOB  would 
have  provided  more  accurate  estimates  of  risk  coefficients  than  OLS-^^.  But 
we  would  not  have  expected  to  find  a significant  difference  in  their  ab- 
solute values.  Kurthciiiiore , one  c.in  easily  see  from  Table  6 that  the. 
difference  is  not  caused  by  just  a fiw  unisiml  estinates  out  of  the  48. 

Every  fractilc  (except  tlie  minimum)  is  lower  for  ROIi  tlian  for  OLS. 

J 
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D . Stronj’.ly  vnryi.nR  pu  ramoters 

The  second  t/pc  of  Individual  case  that  could  be  regarded  as  "unusual" 
Involves  highly  significant  changes  In  the  coefficients.  Table  7 

lists  securities  for  which  estimated  higher-order  polynomial 
coefficients  had  t-ratlos  in  either  the  OLS  or  the  lobust  regressions 
greater  than  3.5  during  the  IbO  wock  s.u.ijile  period.  These  securities  dis- 
played more  significant  non- s tat ionar i t y than  all  others. 

Some  of  these  securities  may  have  entered  this  group  by  chance. 

Whenever  the  extreme  values  of  a sample  distribution  are  ir.ol.ntcd,  sone 
members  of  the  group  will  be  there  simply  from  random  variation.  Wc 
know,  hoviever , from  the  general  tests  for  the  existence  of  non-st at loimi i tv 
reported  in  section  IV-B,  that  some  securities  really  did  have  ctianging 
coof 1 ic lent  s . 

There  are  two  possible  v/ays  to  determine  whether  a jiarticular  securitv 
had  genuine  non-a  tat  ionar  i ty  . First,  if  a polynom  ia.  1 coefficient  is 
significant  In  the  OLS  and  not  the  Robust  regression,  the  probabilitv  is 
greater  that  the  measured  non-stat  ionar  i tv  has  ari.scn  bv  chance.  This 
conclusion  is  based  on  the  fact  that  the  01. S t-ratios  do  not  net  ess.ii  i ly 
conform  to  Student's  distribution  and  arc  not,  thorcfoie,  necessarily 
Subject  to  a known  acceptance-rejection  region.  In  addition,  the 
robust  coefficients  are  more  reliable  because  data  errors,  bad  spots 
on  the  tav>e , or  grossly  abnormal  observations  arc  censored  by  the  robust 
calculation. 

The  second  possibility  for  discriminating  genuine  from  false  cases  1 .s 
a detailed  analysis  of  peculiar  events  for  each  security.  As  we  niontiontd 
early  in  tlio  paper,  them  are  known  influences  on  the  risk  coefficient  and 
these  can  be  examined  directly  for  change.  If,  for  example,  a major  i::suc 
of  new  debt  brought  about  significantly  greater  capital  leverage  during  the 
period  of  observation,  wc  should  find  this  corresponding  to  a positive 
chang.e  In  the  portfolio  risk  coefficient. 
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Tnbic  7 

New  York  Exchnnp.c:  kisLid  Securities  with 
Stionj^ly  Varying  Parameters  1962-(>S 


Kstimacod  Chauc.c  Chaup.e  lit 
in  R) sk, (Robust)  Market  D/E 
'V’VBl.  I'fO  “Vijr”]  Ratio, 

1962-65 


Linear  Significance 

M>  2.0 

oi.s  Ron 


American  CoDunercial  Lines 

American  South  African 
Investment 

CCI  Marquardc 

Crescent  Corp. 

Dan  River  Mills  , 

Dorr-01 iver 
\ 

Florida  Power 

W.R.  Grace 

Madison  Fund 

N.Y.  Central 

Niagara-Mohat/k  Power 

Pacific  Tel  i Tel. 

Peoples  Gas,  Light  S 
Coke 

Proctor  A Gamble 
Richardson-Mcrri 1 1 
Tractor  Supply  Col 
Univ.  Leaf  Tobacco 


-.785 


-.0736 


-1.33 


-1.74  -1.18 


-.862  -.664 


-.  853  -. 


-1 .04 


-1.24 


-.353 


-.872  h. 


-1.90 


-.0935 


-.823 


-.0947 


-.0747 


-2.62 


-.217 


-.036 


-.0019 


.0776 


Zapata  Norncss  Inc 


3.42 


2.04 


-.166 
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33iP.  numbcrt;  in  loblc  7 i;lvc  thi;  tol.il  t;.cim.ncd  ilianf.c  in  llie 
portfolio  ribk  coefflclrnt  from  the  first  date  to  the  last  date  In  the 
sample  period  (July  12,  1962  - July  29,  1965).  Kor  the  4-variable  model, 
which  pormitb  a linear  chanj'.c,  this  is  simply  twice  the  estimated  slope  of 
the  change  with  respect  to  lime. — - lor  the  10  variable  model,  it  is  tVic 
difference  between  the  estimated  time  track  of  the  coefficient  at  the  last  and 
first  dates.  Also  noted  in  the  Table  are  significant  (with  jt[>2.0),  linear 
polynomial  coefficients  (^jj)  f^or  ordinary  least  squares  and  robust  regression. 
This  table  is  our  attempt  to  present  compactly  the  salient  facts  from  the  six 
regressions  and  32  coefficients  estimated  for  these  18  securities. 

The  facts  are  these: 


(1)  lor  12  of  the  18  securities,  there  is  evidence  of  a strong  linear 
change  in  the  risk  cocf  f i c leiil . lliose  securities  have  linear  co- 
efficient absolute  t-ratlos  in  excess  of  2.0  for  both  OLS  and  Ro- 
bust. In  most  of  these  cases,  the  absolute  coefficient  changes  over 
the  sample  period,  estimated  with  the  4-  and  lO-variable  models, 
agree  quite  closely.  Nine  of  12  differ  by  less  than  15  percent. 


(2)  One  security,  Tractor  Supply  Company,  has  significant  01. S Jinc;.r 
coefficients  but  the  corresponding  Robust  coefficients  are  not 
significant  at  all  7—  The  estimated  negative  change  in  the  risk 
coefficient  is  3 times  larger  for  the  4 variable  than  for  the  10 
variable  model  and  the  R is  only  about  .07  in  the  several  re- 
gressions. k'e  tliink  that  tlie  probability  is  high  that  this  security 

,36/ 

had  a spuriously  significant  (>1  S linear  trend. — 


There  arc  five  securities  which  seem  to  have  had  significant  non- 
linear changes  in  risk  that  biought  the  risk  coefficient  hack  close 
to  Its  original  level.  American-South  African,  for  example,  had  a 
trivial  llne.ar  coefficient  but  highly  significant  quadratic,  cubic, 
and  quaitic  terms  wl'ich  iiriply  largo  deviations  of  the  risk  coeffic- 
ient during  tlic  sample  from  Its  beginning  and  ending  values.  (This 
security  is  a curiouslty,  for  its  estimated  mean  risk  coefficient  was 
negative  and  significant  in  all  regressions).  Tor  the  five  securities 
in  tliir.  group,  the  m 1 n imum  difference  between  the  estimated  4-  and  10- 
variable  absolute  cbongos  was  about  3 75  percent.  This  mc.ans,  of  course, 
that  the  4-v.ariablc  model  completely  missed  wli.it  must  li.ivc  been  the 


1/,/ 


A. 
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complex  patli  taken  by  the  coefficient. 

(4)  Tlic  cstin.Hiil  cb.mpes  in  i i '..k  coefficients  are  not  alw.iys  due  to 
capit.il  sttiKlnre  (l,anp,es.  As  ton('.b  support  foi  this  nssettion. 

Table  ^ includes  esticuUed  cbanp.es  in  tbc  debt-e<iuity  ratio  between 
July,  1962  and  July,  1965,  the  beginning  and  ending  sample  points. 

These  arc  market  values  whenever  possible.  Normally,  the  debt 
equity  ratio  is  the  market  value  of  the  bonds  plus  the  book  v.alue 
of  non-inarkc table  debt  divided  by  the  market  value  of  equity  for 
a given  date. 

For  reasons  which  we  cannot  explain,  there  were  no  t-iatios  in  excess  of 
3.5  for  the  intercept  polynomial  coefficients  of  any  sernrity.  This  is 
mild  evidence  against  both  the  Sharpc-l.i  nt  nor  and  Rluck  versions  of  tl.o. 
capital  asset  pricing  model  because  a significant  movement*  in  the  risk 
coefficient  should  be  reflected  in  an  opposite  and  significant  movement 
in  the  intercept  unless  changes  in  the  riskless  rclurn  or  in  the  I'tro-bct.i 
return  just  happened  to  be  offsetting.  The  riskless  rate  or  the  zcio- 
beta  return  cannot  offset  the  risk  coefficient  movements  of  al_^  securities 
because  (1)  some  securities  had  implied  negative  and  scn;.e  had  positive 
changes  and  (2)  some,  of  the  negative  changes  were  associated  with  high 
risk  coefficients  and  vice-versa.—  ^ Coupled  with  this  fact  is  a danger 
that  sampling  variation  was  responsible  for  the  selection  of  some  of  these 
securities  with  highly  varying  parameters.  Thus,  we  do  not  want  to  asset l 
that  a definitive  test  of  any  theory  has  been  provided  by  tbese  last  results. 

In  the  next  section,  portfolios  will  be  used  to  provide  a more  powerful  exami  . 
of  models  which  require  the  intercept  to  move  in  a direction  opposite  to  the 
slope . 

Another  seeming  anom-aly  is  that  none  of  the  84  secuiities  listed  on  the  Air.eiir.’.n 
Stock  Kxchangc  had  t-ratios  as  high  as  3.5  for  any  polynomial  corf  f i cient  r. . 
Howevci  , if  the  probability  of  a security  having  such  a coefficient  is  equal  to 
the  observed  frequency,  18/930  = .00194,  the  probability  is  roughly  10  pcrceirt 
that  none  would  be  observed  in  a random  sample  of  84.  This  probability  would 
not  justify  an  inference  tltat  New  York  and  Amcrictin  Stock  Exchange  listed 
securities  were  truly  different. 
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E.  Disc  rlmt  itJ  1 1 n^;  bcluoon  ag_t,ro|,>u  e and  iiidlvtdual t\onr.  t a 1 1 onar  1 1 y 

We  will  now  discuss  addltlonnl  tests  for  the  existence  of  non- 
i statlonarity  In  the  overall  universe  of  available  data.  In  subsection 

i B,  we  demonstrated  that  our  cxpositional  vehicle,  a sample  of  returns 

on  Individual  assets,  did  Indeed  display  some  kind  of  non- stat Iona r i ty . 

» 

In  the  application  of  our  method  to  on  individual  series,  there  would  be 
. no  need  to  go  beyond  section  B.  Equation  (11)  would  yield  as  much  in- 

i 

! formation  about  the  nature  of  the  non-stat ionarl ty  as  the  single  series 

could  desclose.  In  other  appl ica t 1 ons , however,  as  with  our  data  here, 
the  availability  of  multiple  time  series  might  permit  a finer  degree  of 
inference  about  any  possible  non-stat i onarity  conr.on  to  all  the  series. 

This  could  be  of  considerable  practical  interest.  For  example,  it  might 
not  be  very  interesting  to  a department  store  that  some  customers  were 
departing  in  one  direction  from  a fitted  advertising/sales  model  and 
other  customers  wore  departing  in  the  opposite  direction;  but  it  would  be 
highly  interesting  to  know  if  there  were  a common  change  in  the  coefficients 
of  sales  response  to  certain  kinds  of  advertising. 

Given  many  time  scries,  inference  about  aggregate  statlonarity  and 
d Isor initnation  between  individual  and  aggregate  non-stationarlty  is  indeed 
' possible.  Going  back  to  the  data  we  have  been  using  here,  their  aggre- 

gation would  imply  the  formation  of  portfolios  of  assets  and  the  returns 
on  these  portfolios  would  become  the  object  of  our  aggregate  empirical 
enquiry. 

There  are  actually  two  advantages  to  using  portfolios  rather  than 


individual  securities  in  data  analysis  of  returns.  Tlie  first  ad- 
vantage is  well  known  in  the  asset  literature,  going  back  to  the  original 
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vork  of  Blumc  [1970]  anti  being  refined  by  Black,  Jensen  and  Scholes  [1972] 

•none  others,  Tblf.  Is  to  nverage  out,  by  « cross-sectional  collection,  events 
peenUar  to  Individual  securities.  Ibc  stiuclure  of  any  .IhI  i geneiallug  ,.ro<  c*/. 
can  be  discerned  more  clearly  after  sucli  an  aggregation  because  the  signal 
to  noise  ratio  is  higher.  Of  course,  some  care  mu.st  be  taken  to  assuie  that 
the  portfolios  formed  actually  differ  a priori  in  their  characteristics.  This 

is  the  same  thing  as  saying  that  the  experiment  should  be  designed  to  maXiO.ize 

the  cross-sectional  variation  in  the  parameter  of  interest. 


The  second  advantage  of  portfolio  formation,  which  Is  noted  here  for  the 
first  time  to  our  knowledge  , is  to  smooth  out  eiodcl's  para''ietc.  s 

which  are  unique  to  Individual  firrs.  V.'c  might  find,  for  example,  that  the 
dispersion  of  the  disturhance  tcim  for  a portfolio  l.s  i;iuch  more  stationarv 
than  for  individual  securities.  Ciii  the  other  h.tnd  , cluanges  in  the  runlet 
price  of  risk  or  in  the  riskless  rate  of  interest  would  be  common  to  all  firms 
and  would  thus  affect  portfolios  of  any  sire  to  the  same  degree.  Since  in- 
dividual parameter  non-stationari ly  would  tend  to  wash  out,  the  abilitv  to 
perceive  common  non-stationarity  rust  be  enhanced  by  tlie  portfolio  aggre- 
gation. 


The  number  of  portfolios  foiT'.ed  is  a n.Ttter  of  judgement  . Since  the  t.'^tal 
number  of  securities  is  fixed  (in  our  case  at  930)  the  larger  the  number  of 
securities  included  In  each  portfolio,  the  smaller  nuist  be  the  total  numher 
of  portfolios.  Ideallv,  one  would  lik.c  to  b.ave  a voiv  laige  nu:;l  i of  port- 
folios of  veiy  l.argc  tdzeo.  lie  decided  to  partition  tbe  sample  into  13 
equal-sized  portfolios  of  62  securities  each.  Sccuiitics  were  .a.ssigneJ  to 
portfolios  based  on  the  relative  rank  of  the  average  risk  coefficient,  h 


/ 
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in  model  (11),  from  the  four-var  i.ihl  e toluist  r eg  i es.s  i on  ' cs  t ima  1 1 d over 
the  160  weeks  from  July  5,  1962  t'.reug,h  Julv  22,  1965.  if  . va> 
the  Jth  largest  of  those  930  coefficient  cr.tip.ates,  the  security  w.i«3 


TM’I.K  ’) 


Regress  loiis  on  I'ortfolJos  feu  fv/o  S.ini[>]c  Periods  of 

160  v;oe  r s each 


OLS 

KObUST 

1 Ol.S 

KOBirST 

1962-65 

1965-68 

No.  of 
Var  tables 

1 

Viedian*  value  of  t^ratio  for  average 

risk  coefficient,  b . 

o 

TWO 

22.3 

22. A 

18.3 

19.7 

FOUR 

21.2 

21.5 

18.2 

19.7 

TEN  1 

1 

18.6 

19.2 

12. A 

12.8 

1 

Median*  coefficient  of  determination, 

TWO 

.759 

.678 

.71A 

FOUR 

.756 

.760 

.675 

.712 

TEN 

.757 

.765 

.679 

.713 

i 

* The  average  of  the  7th  and  8th  largest  sample  values  out  of  13 
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Table  10,  a coripaiilou  to  Table  3.  reports  a test  for  tlie  existence 
of  non-statloi'.arlty . It  contains  the  cross-sect  lonal  mean  t-values 
of  the  higher  order  polynomial  coefficients  of  model  (11)  divided  by 
their  estimated  standard  errors.  Only  the  robust  estimates  arc  used 
because  the  sampling  d I.s  t r 1 but  ion  Is  hnoun  to  follov  the  Student  law 
in  this  case,  (see  Section  IV-D)  under  the  null  liypothesls  of  statlo- 
narlty.  Coefficients  t.'lilch  are  significantly  different  from  zero 

provide  evidence  that  the  basic  model  did  not  have  stationary  para- 

43/ 

meters  during  the  sample  period. — The  table  also  reports  the  mean 
coefficients  divided  by  their  standard  errors  as  computed  from  the  cross- 

/A 

sectional  standard  deviations  of  the  coe f f Ic i ent s , (i^  ) . 

Table  10  indicates  that  non-linear  terms  in  the  10-varlable  model 
displayed  the  most  measured  significance  during  the  second  period. 
Although  the  average  estimated  linear  trend  in  the  4-varlahle  model 
was  statistically  significant,  it  was  not  large,  (.041),  and  was 
t either  as  large  in 
linear  coefficients 


absolute  value  nor  as  significant  as  the  non- 

44/ 

in  the  10-varlabIe  model. — 


TABI.E  10 


Tests  for  tlic  Presence  of  non-Stationni  y Paraniolers 
in  15  Portfolios  formed  from  930  NYSE  and  AMEX 
Securities,  1962  - 1968 

(Robust  Estimates) 


Estimated 

No.  of 
Vbls  in 
Mode  1 

1962  - 65 

1965 

- 68 

Pa  rameter 

t/(o^//N) 

c/(o  //n) 
c 

CeneXi:  icl.ejit 

) .(t-  rar.i.g) 

cAo^/>'^l) 

(r  ;T>>  f f i ri />ni: 

Intercept 

Trend, 

4 

3.72* 

6.91* 

.0175 

-1.13 

Slope 

Trend, 

' 4 

-1.56 

-1.85 

2.44* 

3.88* 

Intercept 
Trend,  3^ 

10 

4.88* 

8.34* 

-1.99* 

-3.88* 

Intercept 
Quadratic,  82 

10 

-1.35 

-1.04 

2,94 

5.28 

intercept 
Cubic,  3-j 

10 

-6.77* 

-15.3* 

-8.20* 

-15.2* 

Intercept 

Quartic, 

10 

-5 . 96 

-11.5* 

3.98* 

8.94* 

Slope  ^ 
Trend, 

10 

-.901 

-1.00 

5.17* 

11.0  * 

Slope  ^ 

Quadratic,  ^2 

10 

5.69* 

10.1* 

1.16 

1.35 

Slope  ^ 
Cubic,  bj 

10 

.824 

2,43* 

10.5* 

15.0* 

Slope  ^ 

Quartic, 

10 

- 2.96* 

-3.78* 

-7.06* 

-13.4* 

^Rejection  of  llic  hypothesis  of  stationarity  at  the  95  per  cent 
level  of  significance. 
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There  appears  to  have  been  sljjnlficsnt  changes  In  parameters  during  both 
sub-periods  but  there  seems  to  be  little  connection  between  the  two.  Of 
the  ten  pairs  of  coefficients  computed,  only  two  pairs  are  significant 
with  the  same  sign  and  the  remaining  pairs  have  either  one  non-slgnlflcant 
member  or  are  significant  with  opposite  signs. 

In  comparing  the  portfolio  t-ratlo  results  for  the  first  sub-period  (1962-65) 
with  the  same  computations  for  Individual  securities  (the  next-to-last  column 
of  Table  3}^  the  portfolio  coefficients  are  smaller  In  seven  of  the  eight  cases 
which  were  "signl f leant"  for  Individual  securities.  However,  In  only 
two  cases  did  the  estimated  coefficient  become  "insignificant"  at  the 
95X  level  and  both  of  these  cases  retained  the  same  sign.  As  we  argued 
In  Section  IV-B,  this  test  statistic  will  tend  to  overstate  the  signi- 
ficance of  non-stationarlty  if  there  Is  positive  dependence  cross- 
sectionally  In  the  estimated  coefficients.  The  comparison  of  Tables 
3 and  10  strongly  suggests  that  such  positive  dependence  Is  Indeed  present. 

The  correction  of  this  bias  turns  out  to  be  a rather  difficult  task  be- 
cause there  is  no  a priori  information  about  how  the  covariance  itself  ^ 

might  differ  across  pairs  of  securities  (or  portfolios).  In  order  to 
explain  our  solution  to  this  problem  as  simply  as  possible,  let  us  use 
only  the  simple  two-v.Tr iable  ordinary  le.ist  squares  model  in  devi.ation 
form.  (The  discussion  gcneralizi;s  readily  but  tediously  to  ilic  more 
complex  models  and  tlic  calculations  were  actually  made  for  them) • 

With  the  simple  model  in  deviation  foim,  we  have  for  portfolio  j 


R.  -R.-b^.(R  -R)+  e. 

j,t  j 0,j  ni,t  m j.t 

whore  the  bars  indicate  intertemporal  means.  (Cf . eq.  (1)).  If,  for 

compactness,  we  denote  x = R " R and  the  sample  vectors 
‘ t mt  m ' 

X*  = ]and  rj  = [fj  wlierc  T is  t lie  sompl  c size, 

then  the  least  squares  estimate  of  bjj  j is  given  by 

“ b-  . ♦ (x'x)  ^x'c- 
O.J  O.j  J 
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Now  assiuiiK  tli.it.  fach  1 1 mc  coefficient  (not  its  estimate)  is  (generated 
by  deviations  about  a mc.in  cocfficcnt,  tliat  is 


wticre  is  the  r.ican  ci  oss-sec  i ioua  1 coefficient,  whose  estimation  is 
the  object  of  our  procedure,  and  f. ^ is  a stocliastic  error  for  which  we 
assume  standard  Siihcrical  properties.  Thus,  the  estimated  coefficient 
can  be  considered  the  d( pendent  variable  in  a regression  function  with 
no  explanatory  variables: 


^0  " 


1,...  ,15. 


where 


b.  H , * (x'x)  ^x'e.. 
J 3 J 


The  problem  in  estimating  arises  because  the  di  .stu!  bances  of  this 
model,  the  p's,  arc  not  necessarily  independent.  Hovtcvcr,  one  can 
observe  from  the  model's  structure  that  an  estimate  of  the  complete 
covariance  matrix  of  its  disturbance,  E(yp'),  can  be  obtained  from  a 
combination  of  the  cross-sectional  variance  in  the  estimated  coeffic- 
ients plus  the  estimated  covariance  matrix  of  residuals  from  the  original 
regressions.  In  other  words,  making  the  reasonable  presumption  that, 

C.  and  c.  are  uncorrelated 
J J 

E((kv1j^)  •=  E(Cjf,|^)  (x'x)  ^x'E(£je|^  ) x (x'x)  \ 

with  the  additional  presumption  that  lagged  cross-correlations  of  the 
residuals  arc  uncor rcl a ted , (that  Cov(e.  ,c  .)  = 0 for  i ^ 0) , this 

J,t  K.,t  — 1 

expression  reduces  to 


E(Ujl^)  =•  K(f.C|^)  + (x'x)"’  E(fjC^  ) 


or  in  vector  notation 


34 


L(„p')  “ 1(U,')  ^ i:(tV) 

where  c is  the  (160  x 15)  matrix  of  disiurbances  from  the  15  separate 
recicssions  for  iiulividual  portfolios  (and  tlujs  K(c't)  is  the  con- 
current cross-covariaiicc  matrix  of  d i stui  bantes  ) . 

Since  is  by  assumption  a standard  splicrical  error,  the  of f-diaf;onal 
elements  of  K(66')  arc  zero  and  all  di.ip.onal  clcmcuts  are  the  same  constant, 
denoted  An  estimate  of  is  given,  therefore  by 

(6^)  - (y^)^  I'G  'G1  - l'(x'x)'’  f'e  1 

where  1 is  the  15-element  unit  vector,  e is  the  matiix  of  residnajs 
from  the  15  separate  portfolio  regressions,  and  G r - b^  is  the 
vector  of  deviations  of  coefficiants  ahouL  their  cross-sectional  mean, 


Given  this  estimate  of  the  cross-sectional  variance  of  the  true  co- 
efficient, an  estimate  of  the  covariance  m.atrix  of  disturbances  in  the 
regression  model  above  is 


ft  = £ (Mb' ) I 


+ (x'x)  ' e'e 


where  1 is  the  15  x 15  identity  matrix.  Then  an  estimate  of  the  true 
mean  beta  is  obtaii\cd  from  the  Aitken  estimator  (or  by  generalized 
least  squares ) as , 


b Hd'ft 


and  the  standard  error  of  b is  obtained  as 


s “ I (I’jrh)"' 


Generally  spe. iking,  thi.s  cstir.i.itc  of  can  and  dues  differ  from  the 
simple  cross-sect  lonal  mean  of  the  h^'s,  (wliich  we  denoted  b^  ,ind  used 
to  estimate  the  disLurb.inrcs  G above).  It  is  therefore  natural  to  rc- 
cstimatc  the  u's  using  this  new  and  presumably  more  prcci.se  estimate  of 
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bg.  We  cliu  tliir,  aiul  i li<  n 3tcr;itL(l  the  |.i  or  cdurc  until  tlic  cstim.ilc  b 
chaiif.fd  by  Iom;  thin  .01  I'Oicout.  This  was  accompl  i r.lud  in  .it  most  8 
iterations.  Tht>  rcsultiuj;  t'Stiit.itcs  and  their  t-iatios  using  the 
standard  error  calculated  a.s  above  (but  with  tlic  final  iteration)  are 
given  in  Table  11.  These  t-ratios  are  in  principle  accurate  measures 
of  significance  for  the  .ivci  age  coefficient  after  taking,  into  account 
the  cross-sectional  covariation. 

The  general  impression  of  these  results  on  comparing  them  v/ith  Table 
10,  is  a reduction  in  the  significance  level.  With  one  or  two  exceptions, 
the  signs  remain  unchanp.ed  but  the  absolute  size  of  tbe  estimated  t-ratio 
for  the  higher  order  polynomial  coefficients  is  reduced.  Nevertheless, 
there  are.  still  indications  of  significant  non-st  a t i onai  i t y in  both  periods. 
The  slope  trend  estimates  for  the  1962-65  period  have  actually  increased  in 
significance,  marginally,  in  both  the  6-variablc  and  10-variable  models. 

We  can  be  fairly  certain  that  these  results  are  jivirged  of  all  bias  int  in- 
duced by  cross-sectional  dependence.  Thus  an  inference  of  significantly 
concurrent  non-st.it.  iotutr  iey  In  the  coefficients  of  model  (1)  would  be  very 
reasonable . 
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TAIil.E  11 


Test  s 

fur  the  Presence 
A)  1 k cn-c  St  i lu.i  t or 
t ro  ss-Sec  t i I'lia  ) 

(.f  Non-  St  .It  lo-  T i t y , 
Coi  1 r e i i on  f oi 
Cuv.ri  i at  ion 

ESriH/VlLD  rAUA.’IKTEK 

No.  of  Vbls 
in  tlnde! 

1902  • 09 

1969  (>8  ] 

C - ratios  from  Altken  procedare 

INThKCElT  TRKtJI) 

4 

. 52H 

.481 

Sl.Ol’F.  TREND 

4 

-2.04 

1 .25 

INTERCEl’T  TKKNI) 

10 

1 .4? 

-.89  3 

INTEltCERT  QUADRATIC 

10 

-1.22 

. 

IKTERCEl'T  CU  lUC 

10 

-3.12 

-2.37  I 

INTERCEl’T  QUART  I C 

10 

-2.8/ 

.99  3 1 

SEOl’E  TREED 

10 

-1  .75 

1.92 

SLOPE  QU/VDRATIC 

10 

4.11 

1 .02 

SLOPE  auuc 

10 

-.207 

2.92 

SLOPE  QUARTIC 

10 

* 

-2.26 

-2.17 

Coefficients  and  t ratios  f rt'm  Aitlfii  pi  oc 

ad  Li!  1 1 

MEAN  SLOPE,  b,, 

1.02 

! 

.999  1 

0 

(8.92) 

(12.9)  1 

MEAN  SLOPE, 

10 

1.00 

1 

. 9 -'i  S 1 

0 

(8.93) 

(13.3)  1 

* Rejection  of  t lie  liyfiot  hesi  s of  st.uionarity  at  the  95  percent 
level  of  si  cni  fic.uice  . 


F,  FsJ_rv^j^i^ro)-nt  I nnt.i  to  D.'trct  I’l'fSfiu  i.-  of  Pn  lhn' r r 1 y 1 n{;  Conimpn 
General  Irij*  i'roc  I'-  s 

As  a final  topic  in  tliii:  paper,  vc  wl.'.h  to  prcienl  a sn>jger, 1 1 on  for 
using  the  presence  of  non-stnl lonarl ty  to  examine  the  structure  of  the 
data  generating  process.  Imagine  tfiat  all  of  the  time  series  available  to 
an  Investigator  have  been  generated  by  an  underlying  proccBS  of  identical 
form  but  vith  different  parameters.  Imagine  further  tlmt  the  parameters 
wander  over  time  but  that  tlicir  paths  are  tracked  by  the  methodology  out- 
lined In  the  preceding  sectioiis.  We  assert  that  these  estimated  time  paths 
give  an  opportunity  to  infer  sostething  about  the  generating  process,  if  it 
is  indeed  common.  Intuitively,  the  c los  s- s ec  t i oiia  1 i e 1 .i  t i onsh  i p ami'ng,  the 
tracked  parameters  at  a given  time  shi'nld  have  as  Its  coefficients  the 
values  of  the  factors  In  the  underlying  process  at  that  same  time.  If 
the  form  of  this  c r<'S  s - sec  1 1 ona  1 relation  were  not  s t ,,  t i ona  t y , the  Inves- 
tigator would  he  obliged  to  deny  tbe  validity  of  the  particular  generating 
he  had  hypothrsi zed . 

As  a concrete  example,  we  now  use  the  cross-sectional  sample  of  13 
portfolios  in  a simple  investigation  I'f  the  two-fact  r,  "zero-beta"  gener- 
ating process,  the  process  mentioned  in  footnote  9 in  which  the  coefficient 
equal  to  i 1 h f s model  has  a very  special  stiiic- 

ture,  ns  It  implle.s  that  the  iiUerropl  tatii  jtc  fro'ii  (1)  will  be  g;iven  by 


+ d,  (1-p. 

J.«-  o l.t  ' j.f 


L 


A 


where  E(6q)  “ 0, (and  K(6^)  - the  two-factor  Generating  process  is  valid. 

The  estimated  coefficients  of  (13)  will  be  asynipLot i cal ly  biased 
because  the  explanatory  variable  contains  an  error,  being  an  estimate 
itself.  There  are  two  possible  ways  to  alleviate  this;  (1)  we  could  use 
the  estimated  rcgrcssioit  standard  error  of  to  correct  the  bias  or 

(2)  we  could  use  unstrumcntal  variables  to  obtain  consistent  estimates. 

- . ^5  / . 

Table  9 shows  that  the  t-ratio  for  b^^  is  on  the  order  of  20  — . Since 
the  mean  value  of  b^^  is  on  the  order  of  unity,  the  standard  error  is  on 
the  order  of  .05.  (Over  the  6 models  and  15  portfolios,  the  actual  maximum 
standard  error  was  .0818).  Since  the  cross-sectional  standard  deviation  of 

bp.  is  on  the  order  of  .25,  the  maximum  percentage  decline  in  the  absolute 
value  of  caused  by  the  error  in  bp.  is  around  10.1  percent  and  it  will 
probably  be  less  than  half  this  amount.  The  bias  in  the  intercept  of  (13), 
is  likely  to  be  much  smaller  than  the  bias  in  the  slope.  In  fact,  the 
bias  in  6p  would  be  exactly  zero  it  the  cross-sectional  average  of  bp^  wore 
exactly  unity. 

Results  for  model  (13)  applied  to  the  cross-section  of  15  portolios  are 
given  in  Table  12.  Two  estimates  are  given  for  each  coefficient,  the  first  being 
computed  by  ordinary  least  squares  and  the  second  using  instrumental 
variables.  As  instrumental  variables,  we  used  estimates  of  bpj  for  the 
period  1962-65.  These  were  natural  instruments  since  they  v.’erc  computed 
for  the  same  portfolios  but  ii\  an  earlier  period. 

The  results  indicate  a highly  significant  non-zero  value  for  5p 
which  is  strong  evidence  against  the  two-factor  zero-beta  market  model. 

The  estimates  of  6^  obtained  with  the  instrumental  variables  arc  very 
close  to  those  obtained  with  the  classical  method. 

As  for  the  slope,  supposedly  an  estimate  of  the  moan  zero- 
beta  portfolio  return,  it  is  significantly  negative  in  all  eases.  Notice, 
however,  that  the  effect  of  the  error  in  the  explanatory  variable  shows 
up  in  the  difference  between  ordinary  least  squares  and  instrumental 
variable  estimates.  In  every  case,  the  instrumental  variable  estimate 
is  greater  than  or  equal  to  the  01, S estimate  in  absolute  magnitude.  The 
maximum  difference  is  n.8  percent  f<..r  the  10-variahlc  model  using  the 
OLS  nietliod  to  estimate  bpj . — . This  conforms  very  well  to  the  anticipated 
crrors-in-^ariablcs  bias. 
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Tab 1 c ] 2 

Tcstf;  of  llio  Zt-ro-Bcta  Two-Factor  Model  with  15  I’ortfolioc 

1903-60 


iN'U.MF.FR  OF 
VARIAKI.KS 

2 

/. 

10 



2 

4 

10 

MKTllOi) 

OLS 

p 

OHUST 

6 , INTERCLFT  (Z/ycar) 
(t) 

17.5 

(27.5) 

17.5 

(27.7) 

21  . 1 
(20.1) 

17.7 
(21 .0) 

17.7 

(21.0) 

2 1 . 3 
(22.9) 

6p,  INTKRCtFT 

17.5 

■ T-;.v 

21  .1 

17.7 

17.7 

2 1 . 4 

6 . Sl.OPK  (Z/ycar) 

^ (O 

- 6.0'. 
(2.70) 

-7.1‘' 

( -P.OS) 

- 12.0 
(-3.01) 

-10.1 

(-3.05) 

-10.0 
(-3.2  0 

-15.5 
(-4.01  ) 

6^,  SLOPE 

-6.83 

-7.2'J 

-U.l 

-10.1 

-10.7 

-17.4 

ADJUSTED 

.311 

.355 

.462 

27  2 

.403 

.519 

STUDENT  IZF.D 
RAKCE  OF 
RESIDUALS 

3.55 

3.50 

3.77 

3.77 

3.80 

3.88 

++  Usinf,  as  jnstrvunout.il  variables  the  corresiiondiiij;  estiniatcd  coo  f f i c ; eut  s 
for  1962-65.  (See  text). 


i 
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Tlic  hiphly  sif.ni  f leant  intercept  c.innot  be  cxplnincd  by  errors- 
in-variablcs . In  the  10-vnriabIc  models,  the  bias  would  actually  wore,  in 
the  wrong  direction  since  in  these  eases  tlie  average  value  of  1 - is 

positive  (and  tlius  an  error  adjustment  in  tlio  slope  6^  ; making  it  sipallcr, 
would  move  .6^  even  further  from  zero).  Vor  the  2-  and  U-  variable  m.c.dcls, 
even  if  the  error  in  5^  were  ?5  percent,  (and  it  is  certainly  much  less  th.in 
this),  the  decline  in  the  estincated  intercept  would  be  only  on  the  order  of 
0.1,  leaving  it  essentially  unchanged  and  still  over  20  standard  errors  fiom 

zero . 


We  can  only  conclude  that  the  stochastic  process  generating  asset 
returns  during  1965-68  on  the  New  York  and  American  Exchanges  was  definitely 
not  the  two- factor,  zero-beta  market  mcxiel  with  the  index  we  used.  (Recall 
that  it  was  the  Standard  & Poor's  500  Index).  The  principal  reason  for 


these  results  is  the  exceedingly  bad  perfonaance  of  this  index.  Wtiile 
randcxnly-selected  portfolios  had  average  monthly  returns  in  the  nelgh- 
borhoexi  of  25  percent  per  annum  during  this  sample  period,  the  S&P  500 
had  an  average  monthly  return  of  only  about  9.0  percent.  Since  the  aver- 
age risk  coefficient  of  a randomly-chosen  portfolio  was  near  unity,  (the 

A 

average  for  the  15  portfolios  used  in  this  paper  ranged  between  .995 

and  1.02  over  the  6 regressions  computed),  the  poor  performance  of  the 
S&P  Index  cannot  be  attributed  to  a lower  level  of  "risk"  as  measured  by 
bQ.  This  "risk"  measure  was  not  lower.  Thus,  our  "test"  of  the  two- 
factor  generating  process  is  equivalent  to  the  finding  that  a particular 
wel 1-d iverslf led  portfolio  (the  S&P  Index)  did  not  provide  returns  eexn- 


47/ 

mensurate  with  its  level  of  risk. — 
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It  is  possible  to  use  tlic  entire  time  j>Hths  of  est  Innted  coeffi- 
cients Includlni;  the  h i {;h  or-order  esilr.ntes  of  changes  In  pa  ratncLers 
in  similar  computaC  1 ons . In  fact,  if  it  con  be  presumed  that  I<z , t , 'he 
zero-beta  va li^e,  f ol  1 ows  a random  walV^  with  constant  drift),  then  the  twti 
factor  generatins  process  implies  that  an^^  of  the  higher-order  coefficients 


in  the  equation. 


a..»  6.  + 6.,b.  + £■  ■ , (i^l) 

xj  lO  il  xj  ij 


should  satisfy  the  hypothesis  1,(6^q  ) - 0,  and  ^ 
note— ^ for  proof).  space  considerations,  we  do  not.  report  these 

empirical  results. 


( 
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I 


V.  coNci,i'^:Tor;s 


Vs'e  have  presented  and  explained  a iiiftthodol(n;y  for  using  linear  models 
that  display  non- s t a I I ona  i y pot  ninorers  . Tiie  time  paths  of  slope  and 
constant  terms  were  estii.ated  by  orthigonal  polynomials  of  time.  Noit- 
sta t ionar 1 1 1 cs  in  the  d 1 s t r ibut 1 on  of  the  d 1 sturbances  vere  expurgated 
by  the  technique  of  robust  regression.  The  method  was  applied  to  the 
index  model  of  as.set  returns  (1)  and  the  principal  empirical  results 
were  as  follows; 

1.  There  is  substantial  evidetire  th.at  this  tnodel's  d i ; turhances  are 
non-Caussian  and  there  is  good  evidence  tliat  this  ;s  c.iused,  at 
least  partly,  by  nnn-stat  lon.srity  over  tine  in  the  distribution  of 
disturbances.  The  pethodology  handled  this  problem  adequately  in  the 
sense  that  it  provided  parameter  estimates  '.jlneti  followed  closely 

the  Gaussian  law. 

2.  There  is  strong  evidence  that  the  linear  single  index  model  (i<^ 
non-s t a ti ona ry  parameters.  The  principal  type  of  non-stat i onai i ty 
(or  better  said,  the  most  significant  typo)  seer;s  to  be  common 
across  all  securities.  This  is  consistent  with  the  hypothesis  that 


there  ate  significant  factors  omitted  by  the  simple  model.  There  is 
some  evidence  also  that  non-stationarl ty  is  present  on  the  individual 
security  level.  Since  this  could  not  bo  caused  hy  the  omission  of  a 
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cou,ron  factoi,  it  must  be  due  to  cl....,^es  in  tl.e  c !.a -ac  t c . i s t i c s of  in- 
dividual  firms. 

3.  Tiicie  have  been  episodes  of  unusual  activity  on  the  New  York  and  American 
Stock  Kxchanges.  During  our  sample  period,  1962-65,  the  number  of  out 
liers  (unusu  il  ohsei vat  ions)  clustered  in  the  early  weeks  of  the  record, 
and  there  was  strong  intert cniporal  dependence  in  the  number  of  outliers. 
This  did  not  .seem  to  be  caused  by  variations  over  time  in  the  market  in- 
dex but  rather  by  co-movements  across  securities  in  the  distribution  of 
disturbances . 

4.  The  method  was  able  to  identify  securities  with  particularly  large  char.ge.s 
in  paraa>etcr.s  over  time.  Ttiis  was  true  of  all  parameters.  Securities  that 

have  displayi’d  cxtiomrly  large  rbanp.c.s  in  the  jiortfolio  risk  coefficient 
over  time  were  identified  by  name  as  were  securities  which  bad  unusually 
large  numbers  of  outliers  (and  tlius  displayed  strong  non-stat  i onari t y in 
the  disturbances  or  marked  departure  from  Gaussian  disturbances). 

5.  Refined  investigations  of  non-stationar i ty  and  investigations  of  the 
data  generating  process  were  presented  with  aggregated  data  as  examples 
of  extensions  of  the  basic  method  to  mutlipie  time  series.  Our  method 
can  be  employed  profitably  under  many  circumstances  as  a data  screening 
device,  for  the  Identification  of  individual  cases  with  particular 
non- sta  t ionar  1 r i cs  , or,  it  can  be  used  to  test  hypotheses  which  p.iv- 
port  to  predict  such  changc.k. 


Footnoti'S 


The  theory  of  osset  price  del ei mi n.it ion  was  enunciated  first  by 
Sharpe  [ 1 96''i  1 and  l.i  nt  ner  [1965  ] who  derived  Liarket  cejuilibriuro  con- 
ditions based  on  the  assumption  tliat  all  investors  followed  Markowitz's 
[1952]  advice  to  diversify  their  investments.  In  the  original  Sharpe- 
Lintncr  model,  there  were  only  tv/o  periods,  so  all  parameters  were  sta- 
tionary (trivially).  The  parameter  was  specified  as  aj=(l- 

where  R wa.s  a riskles.s  rate  of  interest.  In  the  later  model  of  black 
F 

[ l973],  however,  there  was  no  riskless  borrowing  and  c<j  bocai.e  relatf?d 
to  the  mean  return  on  a portfolio  whose  return  was  orthogonal  to  K;ri,t- 
Koss  I 197^  has  derived  an  ejui  1 i hr  1 um  model  of  as.set  prices  whiclt  is 
identical  to  the  original  .hliarj  r-Lintner  i:iodel  hut  which  is  based 

on  the  assii  ipt  I on  that  returns  are  generated  hy  (1)  wl  t h^  and  ec'ns  t ant . 
Capital  Asset  Tricing  theory  has  now  come  to  mean  a body  of  ideas  re- 
lated to  risk  aversion,  pottfolio  diversification  and  competitive 
capital  markets.  Furthci  references  (given  below)  will  provide  more 
details  on  the  theory's  liistorical  devclopaient  and  current  cond  i ( i vuir. . 

/ 

See,  for  cx.'.mplc,  Francis  and  Archer,  (1971,  ch.  5],  Sharpe  [ 1970], 
Hossln  11973],  Kama  and  Miller  [1972],  Jensen  [1972],  black[rJ72j 

Some  of  the  events  which  have  been  studied  are  stock  .splltr,  secondary 
offerings,  dividend  changes,  accounting  manipulations,  and  meigeis. 
Brcalcy  [1971]  gives  a reference  list  and  ether  examples. 
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A/ 


4/ 


V 


p/ 


For  an  extensive  analysis  of  this  factor's  influence,  includirif’, 
empirical  testinj^,  sec  Beaver,  Kettler,  and  Scholes  [1970], 
llamada  [1971]  , and  Boncss,  Chen,  and  Jatnsi pi tak  [1974]. 

See  Arrow  [1971,  ch.  3]  . 

There  may  be  other  inflnenres  on  . Kosenherp,  and  McKibhen  [19/i] 

list  over  30  and  present  empirical  evidence  that  at  least  some  of 

them  actually  have  been  influential.  Schaefer,  Brealcy,  Hodptr.  and 

Thomas  [1974]  present  altcrn.Ttive  i"odels  of  the  Bj  as  a stochastic 

process.  Blame  f 1 '> 75 J f I'und  a tendency  for  extreme  values  of  beta 
to  disappear  over  tl:ne. 

Rosenberp  and  McKibben  [1973]  found  six  accounting  measures,  and 
two  historical  measures,  which  seem  to  have  had  significant  in- 
fluence on  the  variance  of  €j. 

The  asset  pricing  theory  of  Merton  [1973]  requires  at  least  one 
additional  explanatory  factor. 

A generalization  (of  the  simple  Sharpe-Lint ucr  model)  recently 
given  prominence  by  Black,  densen,  and  Scholes  [1973]  aitd  Blume 
and  Friend  [1973]  has  the  form 


R.  •=  R (1  -B-)  + 6-  R . + f;  , 

j,t  z,t  ' j J m.t  j,t 


where  R is  knotmi  as  tlic  "zero  beta"  portfolio  return  (due  to  its  linee.r 
z , t 

independence  of  the  market  return,  R^^  ^).  This  model  is  a special  case  of 
(1)  with  stationary  risk  parameter.  It  contains  more  structure,  however. 


because  it  requires  B • •=  6 • "=  constant  and  a. 

J > t J J I 


R , (1 

z ,t 
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Klnf,  [)966]  ';liov«cl  that  liulu:;tty  fnrtors  were  s 1 f 1 can  t in  ex 

plalnin(>,  toir.inon  stock  retuins. 

Evidence  on  tliC  non-Causslan  nature  of  asset  returns  scciss  to  have 
been  p.iven  an  explicit  (and  tliorouph)  treatncnl  first  bv  Mandelbrot 
(1963).  Olhet  evidence  lias  accuniul  .it  ed  over  tiir.e.  .See,  e.p.  Fana 
119G5),  Roll  11970],  Schwarts'  ami  Altman  119/3], 

See  tbe  first  tlireu  references  in  note 

Press  (1967),  Mandelbrot  and  laylor  [1967],  Officer  [1972],  Clark  [197..], 

Barnca  and  Puv;ncs  [1973],  Hsn,  Miller  and  k’ichern  [1973],  Bl-itlberp, 

and  Gonedes  [l97-'<]  c'f<r  altern.itive  explanations  to  the  basic  phenenena 

that  a.sset  leturn.s  are  too  t h 1 ck- t ni  1 cd  to  h.avc  been  generated  by  a 

stationary  Gaussian  process;  i.e.  that  they  have  loo  tnany  extretne 
observations  . 

Data  errors  are  another  source  of  econometric  troubles  which  may 
result  in  seeminply  thick-tailed  residuals  from  an  ordinary  least 
squares  fit  of  mode.l  (I).  In  a recent  paper,  Rosenberp,  and 
Houglcl  1 197'iJ  present  stiikinp  evidence  C'n  the  effect  of  crrois 
and  suggest  tlie  use  of  "a  statistical  method  that  minimizes  the 
effect  of  outliers",  by  discarding  or  truncating  them  or  by  using 
a procedure  "more  'robust'  than  the  quadratic  mi-thods"  (pp.  1308- 
1309).  As  we  shall  sec,  the  methodology  employed  in  this  paper 
has  exactly  this  aim. 

N 

The  characteristic  function  for  a symmetric  stable  distribution 
with  zero  location,  scale  O’,  and  charac tcri s t i c exponent  > is  : 

d ^ (t)  exp  (-  lo  t j^) . 


1 


U.' 

w 

j_9/ 

20/ 


2Ai 

22/ 

23/ 


Fama , Kiitu'r,  .Ifur.t-n  and  Roll  ll969  J 
Kaplan  and  Roll  [l972  .1 
Scholos  Il'J/2j 

Watts  [1973  J.  In  tliis  case,  the  size  of  the  non-zero 
disturbance  seems  to  be  quite  small. 


mean 


A number  of  other  events,  as  well  as  some  of  those  already 

mentioned  above,  ate  suncnaiizcd  in  Bicalcy  [l97l3- 
Cf . Roll  [ 1 97o] 

For  example,  if  the  observed  period  included  two  years  of  weekly 
data  (109  weehs)  , the  time  variable  would  run  -!,  -bl/!>2, 

-50/92 51/52,  I. 

29/  See  Abraniov;itz  and  Stepun  [1909,  ch.  22  3. 

25/  We  wish  to  thank  the  Manaf;ement  Science  Department  of  Vells-Fargo 
Bank  for  making  the  data  available  in  a convenient  form. 

2f]/  The  Thursday  closing  prices  of  each  week  veic  used.  If  Thursday 
was  a holiday,  no  observation  was  used  for  that  week. 

27/  We  began  the  study  in  Pittsburgh  with  an  earlier  version  of  the  tape. 

At  that  point  , we  had  1092  stocks  in  a sample  selected  similarly.  On 
continuinp,  our  work  in  Belgium  on  the  same  tape,  the  records  for  only 
about  100  secuiirics  could  be  read.  Kvidcnlly,  the  tape  had  been 
damaged  in  transport.  Wells-Kargo  kindly  supplied  another  copy. 

On  this  tape,  wc  obtained  tlie  sample  size  of  930. 

It  is  intcicsting,  to  note  that  a number  of  .securities  which  had  been 
read  successfully  on  the  Pittsburgh  system  and  which  the  Belgian 
system  dcsigitatcd  as  faulty  records  on  the  second  tape  were  securities 
with  a very  large  number  of  outliers  computed  by  the  robust  technique. 
As  an  example,  Benguet  had  87  outliers  in  ,'hc  first  160  observations 
in  the  Pilisburgli  runs.  The  IVrlgian  system  rejected  Benguet  as  a 
bad  record.  Wc  tliink  this  is  evidence  that  the  original  tape,  from 
which  both  tapes  were  copied,  actually  had  a bad  spot  in  the  Benguet 


record  that  is  acccplible  on  some  systems.  It  is  finite  likely  that 
Bcnpuct's  fiat  a record  is  free  of  trroi  in  its  original  form,  but 
when  tapes  are  rf>pied  and  reropied,  sliippcd,  anti  nsod  frequently, 
phy.sical  damajp-  tan  occur.  This  j,  upc.is  t r.  tlial  t lie  robust  tecluiique 
is  particularly  v.iluablc  for  l.iip.e  data  banks  tliat  are  distributed 
widely.  Even  if  a piven  macliine  installation  accepts  the  physically- 
datnaped  record  as  valid,  the  technique  may  reject  it  as  an  outlier. 

One  otlier  disquierinp  piece,  of  evidence;  In  the  oidinary  least  squares 
runs  in  Pitlsburph,  tiie  average  standard  deviation  of  regression 
residuals  across  1042  securities  was  10  times  as  high  as  the  average 
for  930  securities  in  tiic  hclgian  runs.  We  later  discovered  that  this 
was  due  to  a single  security!  The  mean  standard  deviations  for  tl>e 
robust  cases  diffeted  only  in  tlie  tliird  decimal  place  in  the  two  runs. 
These  comments  are  very  siir.ilar  in  spirit  to  Uoscuberg's  and  Houglet's 
1)974  ]. 


28/  In  fact,  most  choices  in  empirical  work  are  made  in  such  a sequential 
and  stochasti c nenner;  even  when  the  final  paper  describes  the 
choices  made  as  liaving  been  pre-ordained  bv  logic.  The  actual  pro- 
cedure for  sample  .selcctlcn  and  consol  i tin  t i on  relly  makes  no 
difference  provided  tliat  tiie  results  observed  at  an  intermediate 
step  do  not  influence  the  results  reported  in  the  final  paper. 

The  formula  for  the  asymptotic  standard  deviation  of  a sample 

range  is  obtained  from  known  formulae  for  asymptotic  variances 

and  covariances  of  order  statistics.  See,  c.g.  Cramer  [l946,  p.  369 J 


30/ 


The  logarithm  was  taken  in  order  to  reduce  skcvfncss'  in  the  dependent 
variable. 

y 


j 
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The  rrenu  t v.ilucs  for  tlu“  tot' f I i c ii'iU  n in  Tol’lc  6 ore  as  folli 


No  of 

r i ab  1 e s 

2 

4 

10 

Technique 

OLS 

3.23 

2.94 

2.42 

KOB 

3.39 

3.24 

I 

i 

1 

1 

3 ? / That  is,  estimates  with  hif'.her  t-values. 

33/  It  is  twice  tlie  slope  because  time  lias  been  normalized  to  span 

- 1 to  + 1 . 


34^  In  all  casts,  th'  U variable  ami  10  variable  es  tinaies  ap^reed 
on  the  level  of  significance. 

3 3/ 

The  t-ratios  arc  on  the  order  of  -l.A. 


3f'/  One  othci  piece  of  siip()orting  evidence  for  this  assertion  is  the 

sizeable  difference,  1.2'iS  vei'.iis  .6tl6,  bett.’een  tbe  ri.sk  coefficient 
estimated  (I'ith  the  2,-variable  ncidel)  by  OLS  and  Robust  i cspect  i vely . 
(The  t-ratios  were  5.12  and  3.52  i espect  i ve  1 y ) . This  strongly  sirg./.t'-. 
a bad  and  very  influential  sample  observation. 


All  four  casts  at e ptesent  among  the  18  securities  of  Table  7. 

For  examples  : 

\ 

Security  Direction  of  Risk  Estimated  !!•  an  I'.isk 


Coeff.  Chang.c 

Areer.  Comm.  Lines  -f 

Zapala-Norness  -t 

Florida  Power 

W.K.  Grace  “ 


Coeff.  (lO-VbJ.KOK) 

.801 

1.74 


I 

1 
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Using  f t on  iln;  f oui  - v.i  r i ab  1 c robusL  model  is  also  a somowliat 

arbilraty  clioiro.  We  believe,  ho^;ev(  i , tbal  t lii'  robust  cstinalors 
are  better  than  the  Ol.S  estimators  and  that  some  non-s t a t i ona r i t y 


is  present  in  the  data.  Thus,  only  the  4 and  10  variable  robust 
. estimates  ot  b^j  would  be  good  candidates  for  the  portfolio  assig.n- 
nent  criterion.  The  final  choice  of  the  4-variahlc  rather  than  the 

10-Viiriablc  estimate  was  conplcrely  arbitrary  the  similarity  in 

these  cEtiniates  implies  that  the  choice  did  not  matter. 

39/  The  portfolios  were  formed  with  an  equal  weighting  on  each  security's 
return  in  every  period.  Tliis  is  known  as  an  equally-weighted,  re- 
balanced portfolio  formation  method. 

AO/  During  the  second  sub-sample,  sonic  securities  bad  missing  observa- 
tions. Wlicn  this  occurred,  the  portfolio  for  that  wi’ck  was  con.j’osed 
of  an  equally-weighted  average  of  the  securities  which  d_i_d^  have 
observed  returns.  In  no  case  was  a portfolio  composed  of  tiwer  ih.an 
57  securities  (out  of  62  possible). 


41/  N.B.  The  estimated  risk  ccufficients  for  a portfolio  can  differ  froai 
the.  mean  estimates  for  its  component  securities  v/ith  the  robust  model. 
This  is  because  the  outlier  observation  can  be  and  were  diffeient  wlien 
performing  the  robust  calculation  on  individual  securities  and  cn 
®r.Rtegates  . Thus,  for  example,  the  cross-sectional  mean  it. dividual 
estimates  for  1962»65  were  minimum,  ,20u,  and  maximutu,  1.91  (Cf.T.ibie  b)  , 

7/  The  difference  between  the  two  results  is  highly  significant.  for 

the  individual  securities,  given  that  the  number  of  outlieis  is  com- 

puted from  a t i uc  random  sar.ple,  the  standard  ei i or  of  the  sample  pro- 
portion is  on  the  order  of  .0004.  for  portfolios,  the  associated 
standard  error  is  about  .002.  Thus,  the  difference  in  proportions 
is  at  least  eight  times  its  standard  error. 


n 


AjK>t1ici'  t'  :.t  ini-hiili'd  iti  T.il  If  3 involvfd  Uif  ^ '/'.-r.f'CT  lor^.. 

90  ;.'txen!iif  i-n.p/'-  w ir.  .i  tfst  fur  n -i.-<-'n,rijr-rfnt  nr.u- 

stationaMTv  <ur.  ,i.K  indivi  1,mI  * a.'  sfrie:,).  It  . not 
heix--  N-f.ur"  thc-f  a>f  r.r.ly  Jl.  v a\ie'..  !.■  •)  ;n  tl.e 

cross-sect  i null  distributirn  ui;d  oixler  slatis.tic'  f:sim 
such  a snail  s nri-h'  aa-  r;ot  )-'Ii;i>>Ie.  Tlui  sanjl  inp,  d^stir  i- 
hut  ioji  or'  ' r'lfi'  ■. b at  is.t  ic  • i:  'jnby  •is^nit’pt.ot  ica  1 1 . 

Caur.sian.  Pcside;.,  it  is  quite  likely  that  clianges  in  para- 
metorr.  wliieh  ate  non  lonrur  rent  .icross  scent  it  in*,  would  cancel 
out  in  the  pioc'css  ol  poitK'lif)  ,i  g;;i  ep,.!  t i ni . 'llm.i,  in  .ill 
probability,  they  would  h.rve  been  undetectable  anyway  in  the 
portfolio  d.ibi  , c'.'Ui  if  ll'ey  were  present  in  individual 
sccuiity  iafd. 

,,  Avcrat'C  cross-sectional  tohust  estimates  of  b,  .ind  b for  the 

44  ^ n ) 

15  portfolio;;  duiiiig  the  si  cond  pciiod  were  .3'i0  .ind  -.3P'>, 
respectively.  These  are  the  coefficients  of  the  cubic  and 
quartic  polynomials  of  the  change  in  the  estimated  risk 
coefficient,  denoted  "slope  cubic"  and  "slope  tpiartic"  in 
Tabic  10. 

4^_  For  the  second  sul)-s.ii;;pl  e , hrn.'ever , it  is  only  on  the  ordet  of 
12  for  the  10-vari.ihle  model. 

45,  The  d i f f e 1 c lu’e  s bcli^een  onfinaiy  estiriati  c.  of  6^  and  i ns.I  t ument  1 1 
variable  estimates  suggest  tfiat  the  10-variable  r.iodel  doer,  not 
provide  as  accurate  estimates  of  the  as  the  ?-  and  -'i-v.it  iahle 

models.  However,  the  clostt  agreement  of  tin  dstimite  ■a.'ing 
tlie  robust  method  relative  to  t lie  01. S iiictliod  .ijp.iu  shov;'  ti.'1'’i'.l 


met  hod ' 


supe r loi  i t V . 


Tliis  is  not  oil  th.at  r i ■.  i i.g,  , 'nro  indox'';  whifii  n.iv  .‘■ff.'.'i  tc 
have  very  nimilar  { acpcrt -hui  g.ive  qulf<'  di  1 ; ■ '.'fi.t  tvsul’r.  itr 
rtrvir]  (1).  .‘Vo  J?oll  IldT']  joj  i motv  ccmplef''  inulvair-  of  the 
index'  con;  Of  JU' 'lift  . 


48/  According  to  the  two-factoi  >.  • . >'.i»  ii.e  pr  of 
coefficients  of  (1)  art  g.iven  hy 


02’  ) 


J 


48/  (conLcl)  Define  Lho  v'etor  1’^  rc.nt  aini  n;-,  a:,  c 1 c:i.C’ lU  I lie  or  t hof-.o!!:!  1 

polynomial'-  of  'lil'le  1 i.  v.i  1 lOi  L i .1  ai  t i ii't  1.  ifi.it  . r. , 

P’  H n-  (t)  : r,(t)  ; ...  : P (t)] 

t O 1 n 

wlierc  n is  a nuniber  ;.uf  1 i c i i nt  1 y largo  to  guar.Tritcc  a good 
approximation  - socli  that 


n.  ^ al  P 
J .t  J ' 

■>  ’ J t 

contain  only  trivial  errors  and  ’.•/hor'j  isnd  1;  c^diiern  vprtor:-  cf 

constants  wlvse  estimation  is  the  (ib)cct  of  the  empirical 

fitting, 

. [a„.  . : ...  . , t.!  - [l.„j  . b,.  : ...  , , 

Tfiis  means  that  the  two-factor  rating  rrcicerr.  07')  ran  be  W7':'te^ri 


«•  ^ - O - b!  1'^) 


and  we  note  that  since. 

) 2 for  i = 0 

.(  ’’i  ^ " 0 for  i ? 1 


we  can  intcgiatc  lioth  side.s  of  (13')  to  obt.iin 

a . ° (1  - b„ . ) K wliich  is  eniiivalcnt  to  the  fitted  eqnatioii  (13) 
Oj  Oj  / _ - _ 

in  the  test,  (with  - 0 and  6,  = R ).  This  operation  presnires 
0 1 z 

that  R is  unrelated  to  P.(t)  (i  1)  which  is  satisfied  if 
z , t i 

R flucLn.itC!.  randomly. 

z.t 

To  obtain  a simil.u  test  with  the  full  vector  of  polynomial  co- 
efficients, wc  defin.'  the  in  t e gr.i  t ed  outer  pioduct  m.ntrix  of  P^ 
with  itself  as 
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